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A fully geometrical treatment of general relativistic magnetohydrodynamics (GRMHD) is devel- 
oped under the hypotheses of perfect conductivity, stationarity and axisymmetry. The spacetime 
is not assumed to be circular, which allows for greater generality than the Kerr-type spacetimes 
usually considered in GRMHD. Expressing the electromagnetic field tensor solely in terms of three 
scalar fields related to the spacetime symmetries, we generalize previously obtained results in various 
directions. In particular, we present the first relativistic version of the Soloviev transfield equation, 
subcases of which lead to fully covariant versions of the Grad-Shafranov equation and of the Stokes 
equation in the hydrodynamical limit. We have also derived, as another subcase of the relativistic 
Soloviev equation, the equation governing magnetohydrodynamical equilibria with purely toroidal 
magnetic fields in stationary and axisymmetric spacetimes. 

PACS numbers: 04.20.-q, 04.40.Dg, 04.40.Nr, 52.30.Cv, 95.30.Qd 



I. INTRODUCTION 

General relativistic magnetohydrodynamics 
(GRMHD) is a rapidly developing field of modern 
astrophysics [J-O, driven by numerous observations of 
accretion disks around black holes [JL jets in active 
galactic nuclei or microquasars [^jMSi gamma ray 
bursts, hypernovae, pulsars [l| and strongly magnetized 
neutrons stars (magnetars). In a first approximation, all 
these systems are stationary and axisymmetric. While 
GRMHD had been formulated by Lichnerowicz in 1967 
0, its development for stationary and axisymmetric 
spacetimes originates in the work of Bekenstein and 
Oron (1978) (1 (hereafter BO) and Carter (1979) 0. 
In particular, BO have established two conservation 
laws associated with the spacetime symmetries, the 
first of them being a generalization of the Bernoulli 
theorem to the case of a magnetized fluid. Another 
important step has been the GR generalization of the 
famous Grad-Shafranov equation to the Schwarzschild 
spacetime by Mobarry and Lovelace (1986) [l3| and to 
the Kerr spacetime by Nitta, Takahashi and Tomimatsu 
(1991) [nl and Beskin and Pariev (1993) [H. The 
extension of the Grad-Shafranov equation to the most 
general stationary and axisymmetric spacetimes has 
been performed by loka and Sasaki (2003) [l^l, most 
general meaning without the assumption of circularity 
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(also called orthogonal transitivit''^ , which holds for the 
Kerr spacetime. 

All the studies mentioned above either (i) involve 
coordinate-dependent quantities or (ii) introduce some 
extra-structure in spacetime, such as foliations by 2- 
surfaces, in addition to the canonical structures induced 
by the two spacetime symmetries (stationarity and ax- 
isymmetry). For instance, two of the fundamental quan- 
tities introduced by BO are defined in terms of the 
components Fa^ and of the electromagnetic ten- 
sor and the fluid 4- velocity by uj := — foiZ-Fai and 
C := F^i / {y/ — gnu^) . From these expressions, it is not 
obvious that these quantities are actually coordinate- 
independent. Another example, related to the feature 
(ii) mentioned above, is the (2+l)-|-l decomposition [l^ 
used by loka and Sasaki [l^ [l5| in their study of noncir- 
cular spacetimes. 

In this article, we undertake a systematic study of sta- 
tionary and axisymmetric GRMHD relying solely on the 
spacetime structure induced by the spacetime symme- 
tries. To this aim, we make an extensive use of Cartan's 
exterior calculus, relying on the nature of the electromag- 
netic field as a 2-form and the well known formulation of 
Maxwell's equations by means of the exterior derivative 
operator. We also employ the possibly less well known 
formulation of hydrodynamics in terms of the fluid vor- 
ticity 2-form, originating in the works of Synge fl6j and 
Lichnerowicz 17| . This enables us to formulate GRMHD 
entirely in terms of exterior forms. Such an approach 
is not only elegant and fully covariant, but also makes 
easier some calculations which turn to be tedious in the 
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component approach. We pay attention to keeping hy- 
potheses to a strict minimum, which allows us to present 
the results in their most general form, including noncircu- 
lar spacetimes, and to encompass some special cases that 
had not been considered before, in particular those corre- 
sponding to a pure rotational fluid motion (no meridional 
circulation) or to a purely toroidal magnetic field. 

The plan of the article is as follows. In Sec. we es- 
tablish the most general form of a stationary and axisym- 
metric electromagnetic field, independently of any MHD 
context. In Sec. lIIIi we introduce the concept of a perfect 
conductor and in Sec. lIVI that of a perfect fluid, leading to 
the MHD-Euler equation. We also derive two Bernoulli- 
like conservation laws in that section. Section |V] is de- 
voted to the integration of the MHD-Euler equation by its 
reduction to the master transfield equation, a relativistic 
generalization of the Soloviev transfield equation. Vari- 
ous subcases of that equation are examined in Sec. I VII 
making the link with preceding results in the literature. 
Finally, Sec. I VII] provides a summary and concluding re- 
marks. 



II. STATIONARY AND AXISYMMETRIC 
ELECTROMAGNETIC FIELDS 

A. Framework and notations 

We consider a spacetime i.e. a four- 

dimensional real manifold ^ endowed with a Lorcntzian 

metric g of signature ( V ++). We assume that 

is orientable (see Appendix |B]) , so that we have at our 
disposal the Levi-Civita tensor e (also called volume el- 
ement) associated with the metric g [cf. Eq. (|B1I) ]. Let 
V be the covariant derivative associated with g: Vg = 
and Ve = 0. 

We shall mostly use an index-free notation, denoting 
vectors and tensors on M by boldface symbols. Given a 
vector we denote by v_ the linear form associated to v 
by the metric tensor, i.e. the linear form defined by 

v:=g{v,). (2.1) 

Besides, given a linear form a;, we denote by u) the vector 
associated to by the metric tensor: 

u;=:g(a;,.). (2.2) 

In a given basis (e^), where the components of g, v and 
u> are t/a^, and Wq respectively, the components of 
and ijj are Va = QayiV^ and = g°''^uj^. 

Given a vector v and a tensor T of type (0, n) {n> 1), 
i.e. a n- linear form (a linear form for n = 1, a bilinear 
form for n = 2, etc.), we denote hy v ■ T (resp. T • v) 
the (n— l)-linear form obtained by setting the first (resp. 
last) argument of T to v: 

t;-T;=r(i;, .,...,.) (2.3a) 
T-v:=Ti.,...,.,v). (2.3b) 



Thanks to the above conventions, we may write the scalar 
product of two vectors u and v as 

g{u, v) = u- V ~ u ■ v_. (2-4) 

We denote by V- the covariant divergence, with contrac- 
tion taken on the last index. For instance, for a tensor 
field T of type (2, 0), V • T is the vector field defined by 

V-T:= V,,T"^e„, (2.5) 

where (cq,) is the vector basis with respect to which the 
components V^T"^ of VT are taken. Note that the 
convention for the divergence does not follow the rule for 
the contraction with a vector: in (j2.3ap the contraction 
is performed on the first index. 



B. Stationary and axisymmetric spacetimes 

We assume that the spacetime (^, g) possesses two 
symmetries: (i) stationarity: there exists a group ac- 
tion of (R, -|-) on ^ whose orbits are timclike curves and 
which leaves g invariant; (ii) axisymmetry: there exists 
a group action of SO (2) on ^ whose fixed points form 
a 2-dimensional submanifold A C ^ and which leaves g 
invariant (see e.g. Ref. [l^ for an extended discussion). 
A is called the rotation axis. To each parametrization of 
the one-dimensional Lie groups (M, and S0(2), there 
corresponds a parametrization of the action orbits; the 
corresponding tangent vector fields, called the generators 
of the symmetry group, are denoted ^ for (R, -I-) and x 
for S0(2). The invariance of the metric under the actions 
of (M, 4-) and SO (2) is translated by the vanishing of the 
Lie derivative of g along each generator: 

C^g^O and C^g = 0. (2.6) 

The definition of the Lie derivative is recalled in Ap- 
pendix 1^ Thanks to the expression (|A8|) and to the 
identity V^^q^ = 0, Eqs. (|2.6p are equivalent to the so- 
called Killing equations: 

V„?;3 + V^Ca = and V«x/? + V^Xa = 0. (2.7) 

The group generators ^ and x are then called Killing 
vectors. For a given group action, a Killing vector is 
defined up to a constant factor, corresponding to the 
change of parametrization of the group. Regarding the 
S0(2) action, we can specify uniquely the Killing vec- 
tor X by demanding that it corresponds to the standard 
parametrization of the group S0(2), i.e. by selecting the 
parameter as being the rotation angle ip € [0, 27r[. For 
the (M, -I-) action, there is a priori no natural scaling of 
the parameter < e K. But if the spacetime is asymptot- 
ically fiat, we may fix the scaling by demanding that ^ 
has the standard normalization at infinity: 

^•^^-1. (2.8) 
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For a spacetime that is both stationary and axisym- 
metric, Carter [l^ has shown that no generality is lost by 
considering that the stationary and axisymmetric actions 
commute. In other words, the spacetime g) is invari- 
ant under the action of the Abelian group (M, +) x S0(2), 
and not only under the actions of (R, +) and SO (2) sep- 
arately. It is equivalent to say that the Killing vectors 
commute: 



(2.9) 



Thanks to the property p.9p , one may introduce coordi- 
nates (x") = {t,x^,x'^,ip) on ^ such that 



d_ 

dt 



and 



d_ 

dip 



(2.10) 



Such coordinates arc called adapted to the spacetime sym- 
metries. Within them, the metric components are func- 
tions of {x^jx"^) only: 



9aP = gal3{x^,x'^)- 



(2.11) 



Adapted coordinate systems are by no means unique: any 
change of the type 

t' =t + Fo{x^,x^) (2.12a) 
x'^=Fi{x\x^) (2.12b) 
x'^^F2{x\x^) (2.12c) 

^' ^^ + F3{x\x^), (2.12d) 

where Fa{x^,x'^) are well behaved functions, leads to co- 
ordinates that are still adapted to the spacetime symme- 
tries. 

Usin g th e same notation as Carter in his Lcs Houchcs 
lecture [2^, we introduce the following scalar fields: 



V 

w 

X 



i-x 

x-x 

:= - det 



x-i 



ix 

X ■ X 



(2.13) 
(2.14) 
(2.15) 

VX + W"^. (2.16) 



Since ^ is assumed to be timelike, we have V > Q. Be- 
sides, since % is spacelike, X > 0, except on the rotation 
axis A where X ~ Q. Consequently, cr > except on A, 
where cr = 0. For the Minkowski spacetime and using 
adapted coordinates (i, r, 0, Lp) of spherical type. 



Mink.: 



F = l, W = Q, X 



7'^ sin^ ( 



(2.17) 



We shall also need the Newtonian values of these fields 
to take non-relativistic limits. In standard isotropic co- 
ordinates, 



Newt. 



1 + 2$ 
(1 



W = Q 



(7 = r sm 



grav 1 

2$grav)?'^ sin^ e 

2 , 



(2.18) 



where $giav is the Newtonian gravitational potential, 
which obeys |$grav| ^ 1- Note that throughout the arti- 
cle, we are using units such that c = 1. 



C. Orthogonal decomposition of the tangent spaces 
and circular spacetimes 

The properties of stationarity and axisymmetry de- 
fine privileged 2-surfaces in spacetime: the surfaces 
of transitivity of the group action (R, x S0(2). They 
are spanned by coordinates {t, (p) of the type (|2.10p and 
the Killing vectors (^, %) are everywhere tangent to them. 
Except on A, (^, %) constitutes a vector basis of the 2- 
planc n tangent to S^: 



n = Span(|,x) 



(2.19) 



The metric induced by g in the 2-plane 11 being non- 
degenerate (n is a timelike plane), the tangent space 
,%:{^) at any point x G ^ can be orthogonally de- 
composed as the direct sum 



n® n 



(2.20) 



where 11-'- is the (spacelikc) 2-plane orthogonal to 11. A 
vector V £ ,%:{^) is said to be toroidal iff e 11 with a 
non- vanishing component along x a-nd poloidal or merid- 
ional iff V £ n-*-. 

A question that naturally arises is whether the decom- 
position (j2.20p is integrable, i.e. whether there exists a 
family of 2-surfaces such that at every point H-*- is tan- 
gent to one of these surfaces, in the same way as the 11 
planes are everywhere tangent to the surfaces of transi- 
tivity .y. The spacetimes for which this property holds 
are called orthogonally transitive or circular |20l . |21| . Ac- 
cording to the Frobenius theorem (see e.g. Appendix B 
of Ref. [12] or Sec. 7.2 of Ref. [1^), the necessary and 
sufficient conditions for (./#, g) to be circular are 



and 



0, 



(2.21) 



where 'to^ and are the two twist scalars defined by 

'^i ■■= *(i A X A d|) = e^^'-P-^^x.^p^a (2.22a) 
'^x Ai^XAdx) = e^'^p'^^^x.^px.- (2.22b) 

In these equations, d is the exterior derivative, A the 
exterior product and 7k- the Hodge star; all these operators 
are defined in Appendix |B] 

If (^, g) is circular, one may choose the adapted coor- 
dinates (t, x^, x^, (/j) so that (x^jX^) span the 2-surfaces 
orthogonal to the surfaces of transitivity. This leads to 
the following simplifications in the components of the 
metric tensor: 

.901 = .902 = .931 = .932 = (circular spacetime). 

(2.23) 

Examples of circular spacetimes are the Kerr-Newman 
spacetime (cf. Appendix |D]) and the spacetime generated 
by a rotating fluid star with a purely poloidal magnetic 
field [13, [li] or a purely toroidal one [H [l^- In this 
article, we do not restrict ourselves to the circular case. 
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D. Stationary and axisymmetric electromagnetic 
field 



We consider an electromagnetic field in described 
by the electromagnetic 2-form F, which obeys Maxwell's 
equations: 



dF 



(2.24) 
(2.25) 



where d is the exterior derivative (cf. Appendix [BJ , -kF 
is the 2-form Hodge-dual of F [cf. Eq. (|B5)) ]: 



(2.26) 



★ji is the 3-form Hodge-dual of the 1-form j associated 
with the electric 4-current j [cf. Eq. (|B7p ] : 



*J := e(j', ■) 



(2.27) 



and /xo is the vacuum permeability. 

We assume that the electromagnetic field is both sta- 
tionary and axisymmetric. This is expressed by the van- 
ishing of the Lie derivatives of the electromagnetic tensor 
along the symmetry generators ^ and X'- 

C^F = and C^F^O. (2.28) 



Now, thanks to the Cartan identity (|B2ip and the 
Maxwell equation ((2:241) £ g F = ^ • dF + d(^ • F) = 
d(^ • F). Hence Eqs. (|2.28p are equivalent to 



d(| • F) = and d{x ■ F) 



0. 



Invoking the Poincare lemma, we conclude that there ex- 
ist (at least locally) two scalar fields $ and ^E* such that 



I • F = -d$, 



(2.29) 
(2.30) 



$ and 5" are defined up to some additive constant and 
the minus sign is chosen for later convenience. 

One very often introduces an electromagnetic 4- 
potential, i.e. a 1-form A such that F = dA. Thanks to 
the identity ddA = [cf. Eq. (|BT7)) ]. the MaxweU equa- 
tion (|2.24[) is then automatically satisfied. As shown in 
Appendix [Cl one may use the gauge freedom on A to set 



^ = A- $ = At and 'i' = A ■ x = 



(2.31) 



where the equalities with At and A^ rely on an adapted 
coordinate system (t, x^, x^, </?) [cf. Eq. (j2.10p ]. In this 
article, instead of A we will use the gauge-independent 
quantities $ and 

From Eq. ([2211), $ = 4 • d* = -F(|,^) = and 
C^<^ = X - = -F{^,X)- Similarly, from Eq. ((OOl) . 
C^^ = X d^ = -F{x,x) = and = ^ • d* = 

F{^7 x)- Hence we have 

d[F(^, X)] = d[| • d*] - £^ d* = (F • x) = 0, 



from which wc conclude that F(^, x) is constant over 
We assume that this constant is zero. In particular this 
is the case if F vanishes at some place (e.g. at spatial 
infinity) : 



F{tx)=0- 



(2.32) 



A consequence of the above property is that the poten- 
tials $ and ^E* obey both spacetime symmetries: 

£^$ = £^$ = and £^ * = Z:^ * = 0. (2.33) 

Apart from F{^,x), the only non trivial scalar that 
one can form from F, ^ and x is 



I:^*F{tx)- 



(2.34) 



We may then assert that the most general form of a sta- 
tionary and axisymmetric electromagnetic field is 

F^d'fAC+d^ Ax* + -e{tX,;-) (2-35) 

cr 



*F = e(V$, e, •) + ^Vl-, X*, -,■)-- I A X, (2.36) 

a — — 

where the 1-forms (^*, x*) constitute the dual basis of the 
vector basis (4, x) of the plane H defined by Eq. p.l9p 
and vanish on H's orthogonal complement: 

r-|=l, r-X = 0, X*-^ = 0, X*-X=l, (2.37) 



VweH^, ^* - v = x* -v = 0. (2.38) 



Conditions p.37p - p.38|) define (^*,X*) uniquely. Indeed 
it is easy to see that, in terms of the scalars defined by 

e = -{~X^ + Wx) and X* = -iW^ + Vx)- 
(7 — — a — — 

(2.39) 

In terms of coordinates (t,x^,x'^,(p) adapted to the space- 
time symmetries, wc may express the 1-forms ^* and x* 
as 

e=dt+- i-X^a + Wxa) dx'^ (2.40a) 
a 

X* =d^+i(I^ea + Vxa)dx^ (2.40b) 
a 

where the index a ranges from 1 to 2. In particular, for 
circul ar spacetimes, £,a = 9a0 = and Xa = 9a3 = [cf- 
(I22S1)], so that 

^* = dt and x* = d(p (circular spacetime). (2.41) 

To demonstrate (|2.35p let us consider the 2-form 

G := F - d$ A ^* - d* A X*- 
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It satisfies 

G(i,.) = F(i,.)-(|_^)r 

-d* 

+ (x^)d* = 0. 



Similarly, G{x, ■) = 0. Hence the 2-forni G vanishes on 
the plane H, i.e. the non-trivial action of G takes place 
in the plane H^. Another 2-form that shares the same 
properties as G is 



H 



1 



Since the vector space of 2-forms in the 2-plane 11-'- is of 
dimension 1, G\jj± and i?ln-L must be colinear. Since 
the H\yi± is not vanishing, we conclude that there must 
exist some coefficient / such that G\yi± =I H\jj±. Since 
both 2-forms vanish on 11, we may extend the equality to 
G and H, thanks to the property (|2.20p : 

G = IH. 



This proves that F takes the form (|2.35p . Using the 
properties (|B12[) . it is then immediate to show that the 
Hodge dual of F is given by (|2.36p . On this form, we 
verifj0 that *F{^,x) = I: i-©- that the proportional- 
ity coefficient / introduced above is indeed the quantity 
defined by Eq. (|2.34p . This completes the proof of the 
decomposition (|2.35l) of F. 

Equation (|2.35p shows that a stationary and axisym- 
metric electromagnetic field is entirely described by three 
scalar fields: and /. A concrete example is provided 
by the Kerr-Newman electromagnetic field presented in 
Appendix |D] The component expression of (|2.35|) with 
respect to an adapted coordinate system is given in Ap- 
pendix [e] 



E. Maxwell equations 



The first Maxwell equation, Eq. (|2.24p . is automat- 
ically satisfied by the form (|2.35p of F, whatever the 
values of $, and /. Indeed, since dd$ = dd^* = 
[Eq. (|B17P ]. we have, using the Leibniz rule (|B18p with 
p=l, 

dF = -d$Ad^*-d1'Adx*+d[/CT-ie(|,X, .,■)] (2-42) 

and each of the three terms in the right-hand side van- 
ishes. Regarding the first term, we have, via the Cartan 
identity. 



0. 



2 using the fact that e( V$, ^* , x) = 0, thanks to relation 1 12.391 1 
which induces that the vector ^* is a linear combination of ^ and 



Similarly, x ■ d^* = 0. Hence the 2-form d^* vanishes on 
H. The same thing holds for the 1-form d$, by virtue of 
Eq. (|2.33p . Consequently, the 3-form d$ A d|* vanishes 
on H and acts only in H-'-. Since dimH-'- = 2, the 3-form 
d$ A d^* necessarily vanishes on H-'-. We thus conclude 
that d$ A d^* = in all space. The same property holds 
for the 3- form d^* A dx* ■ Finally, regarding the third 
term in (|2.42p . let us take its Hodge dual and write, using 

*d[Ia-'eit X, ., .)] = - *d * {la-' | A x)- 

Now, the operator -^d-k is the codifferential and can be 
expressed as the divergence taken with the V connection: 
id -k {Ia~' I A x) = V • {Icj-^ i ^ ^o'^ it is '^asy to 
see that V- (/cr^^ ^ A x) = la^' x] = by virtue of 
Eq. (j2.9p . Hence ★d[/o-~^e(^, x, ., .)] — 0, which implies 
d[/o-"ie(|,X, •, •)] = 0- We conclude that ((^^ reduces 
to dF = 0, i.e. the first Maxwell equation dUU). 

The second Maxwell equation, Eq. (|2.25p . gives the 
electric 4-current j. Let us first fix the first two argu- 
ments of each 3-form appearing in Eq. (|2.25p to (^,x)- 

/^oe(j,^,X,.) = d*F(ex,.)- (2.43) 
Now, by means of the Cartan identity, 

d*F(^, ., .) = *F -d[*F(^, .)] = -d[^(^, .)]. 



Hence 

d*F(C,X,-) = -X-d[*F(|,.)] 

= -{C^ [^(t.)]-d[^(|,x)]}=d/. 



Therefore Eq. ()2.43p becomes 

A*oe(i,^,X, •) = dl. 



(2.44) 



We conclude that if the 4-current has some poloidal part, 
i.e. if j ^ H, then necessarily / 7^ 0. 

Taking the Hodge dual of (|2.44p and applying the re- 
sulting 3-form to the couple x) yields 



3 = {i* • j)4 + (X* ■3)X 



1 



■e(|,X,V/,.). (2.45a) 



There remains to evaluate ^* j and x* j- To this aim, let 
us consider the Maxwell equation (|2.25p in its dual form 
V- = Hoj. Substituting (|2.35p for in it, expanding 
and making use of (|2.39p and p.7p results in 



X. 



w 



A^oX 



fW V 
-Vu — V^$ + — V^* 

\ (7 O 



(2.45b) 



(2.45c) 
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where the twist scalars and are defined by (|2.22|) . 

At the Newtonian hmit, the electric 4-current j can 
be decomposed into the charge density pe and the elec- 
tric 3-current J, both measured by the stationary ob- 
server, according to j = p^i, + J and ^ ■ J = 0. 
From Eq. (|2.45a[) . wc get Pc ~ ^* ■ j and J = {x* ■ 
j)x — A*of •)• Choosing spherical coordi- 

nates {t, r, 6, (f ) and using Eq. (|2.18p as well as Eq. (|E1[) 



with 



9 



' sinO to express e(^,x, V/, .), we obtain 



the Newtonian limit of Eqs. p.45|) as 



Newt. 



1 



sin 9 



(2.46a) 



Newt.: IIqJ 



r smf 



-dole 



(r) 



dri e 



(0) 



(2.46b) 

where {e(^r), ^{e), ^{ip)) is the standard orthonormal basis 
associated with spherical coordinates: e^^) :~ dr, 6(0) := 
r~^de and Bf^^-j := {r sm6)^^dip, and A* is the second- 
order differential operator defined by 



sin 9 , 



1 



(2.47) 



III. PERFECT CONDUCTOR 
A. Definition and first properties 

From now on, we assume that a part V C ./# of space- 
time is occupied by a perfect conductor. By this, wc mean 
that D is covered by a congruence of timclikc worldlincfH 
such that the observers associated with each worldlinc 
measure a vanishing electric field. This expresses the in- 
finite conductivity condition via Ohm's law. Let us recall 
that the electric field f-form e and the magnetic field vec- 
tor b measured by an observer of 4-velocity u are given 
in terms of F by 



e ^ F ■ u and b ~ u ■ iF". 



(3.1) 



Equivalent ly, F is entirely expressible in terms of e, 6 
and It as 



F = M A e + e(M, 6, ., .) 
-kF = —u Ab + e{u, e, ., .). 



(3.2a) 
(3.2b) 



The perfect conductor condition is e = 0. From (|3.ip . 
this is equivalent to 



F • M 0. 



(3.3) 



^ later on, we will specify these worldlines to be those of a perfect 
fluid, but this not necessary for the present discussion 



The electromagnetic field then reduces to 

F = e(M,5, .,.) (3.4a) 
*F = -M A b. (3.4b) 

Let us decompose the 4-vclocity u orthogonally with 
respect to the 2-plane 11, thereby introducing the scalars 
A and ft and the vector w : 



(3.5) 



The 4-velocity normalization relation jx -it = — 1 is equiv- 
alent to the following relation between A, H. and w [cf. 
Eq. (j2l^ - (j2l6| ]: 



A 



I + w ■ w 



V - 2wn - xn^ 



(3.6) 



For circular spacetimes and in coordinates {t,x^,x'^,ip) 
adapted to the spacetime symmetries, one has w'^ = = 
[Eq. (Hill below], so that d(p/dt = u^/u° = n, show- 
ing that Q is the angular velocity of the fluid about the 
rotation axis. On the other side, we shall call w the 
meridional velocity. We shall say that the fluid is in pure 
rotational motion iS w = 0; the 4-vclocity u is then a 
linear combination of the two Killing vectors. 

Cartan's identity (jB21[) . along with Maxwell equation 
(|2.24p and the perfect conductor condition p.3p gives 
CuF = u-dF + d{u ■ F) =0 + 0, i.e. 



C^F = 0. 



(3.7) 



This result is independent of the hypotheses of station- 
arity and axisymmetry and is the geometrical expression 
of Alfven's theorem about magnetic flux freezing. 

From the very definition of $ and 5* [Eqs. (|2.29p - 
(jOD)) ]. wc have $ = m • d$ = -F(|, u) and £^ * = 
-~F{x,u). The perfect conductor condition (|3.3p gives 
then immediately 



/:„, $ = 



and 



(3.8) 



Hence the potentials <f> and are preserved along the 
fluid lines. The expansion (j3.5l) of u and the symmetry 
properties (|2.33p show that p.8p is actually equivalent to 



C,„ $ = 



and 



£,„ ^ = 0. 



(3.9) 



Let us express the perfect conductor condition (|3.3p by 
replacing F by its expression (|2.35p : wc get 

{C ■ M)d$ - (d$^)|* + (x* • u)d^ 

^(d^^)x* + ^e(C,X,.,M) = 0, 

-E(C,X-t",-) 

where use has been made of p.9p . Since A ^ (otherwise 
u would be null), we obtain 

d$ = -r!d«' + -^e(^,X,t«,.)- (3-10) 
crA 
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B. Conservation of baryon number and stream 
function 



If n denotes the baryon number density in the fluid 
frame, the law of baryon number conservation is V • 
[nu) = 0, or equivalently, thanks to the decomposition 
p.5|) with ^ and x Killing vectors, 



V- (nw) = 0. 



(3.11) 



Thanks to the identities (|B20p and (|B7|) , we may rewrite 
the above property as 



d(n -kw) = 0. 



(3.12) 



From the Poincare lemma (cf. Appendix[B]), we conclude 
that there exists a 2-form H such that 



dH = n i^w. 



(3.13) 



The above relation is analogous to Maxwell equation 
(|2.25|) . via the identifications -kF o H and fioj o nw. 
Consequently, the same reasoning that led to Eq. (j2.44p 
results in 



(3.14) 



where the scalar field / is related to H by the analogue 
of Eq. p.34p : / := H{$, %). We also have the analogue 
of Eq. (|2.45ap , with ^* - w = and x* ■ if = in addition, 
since w e 11^: 



w = — ^e(^,X, V/, 



(3.15) 



This relation shows that the fluid meridional velocity w 
is entirely described by the scalar field /; / is called 
the stream function (or Stokes stream function). From 
Eq. p.l4p . we have immediately $-df = and x-d/ = 0, 
which shows that / obeys the two spacetime symmetries. 
Moreover, a direct consequence of Eq. (|3.15l) is w df ~ 0. 
In view of Eq. (|3.5p , this yields 



f = 0, 



(3.16) 



i.e. / is conserved along the fluid lines. 

The Newtonian limit of Eq. p.isp is easily taken via 
Eqs. (prS)) and ([ET]) : 



Newt.: 



1 



W = 



nr sin 6 



-def e(r) 
r ^ ' 



drfe^e) 



(3.17) 



where the notation is the same as in Eq. (|2.46p . 

Taking (|3.14p into account, the perfect conductor rela- 
tion ()3.10|) becomes 



d/. 



(3.18) 



Thanks to Eq. p.lSp . the condition ii; • d<i> = 
[Eq. p.9p ] is equivalent to 

e(l,X,V/,V$) = 0. 



This relation is satisfied if, and only if, the 1-forms df 
and d$ are linearly dependent, i.e. if there exist some 
scalar fields a and (3 not simultaneously vanishing such 
that 



ad$ + /3d/ = 0. 



(3.19) 



Similarly the condition iD-d^' = [Eq. (|3.9p ] leads to the 
existence of two scalar fields a' and /?' not simultaneously 
vanishing such that 



a' d* + /3'd/ = 0. 



C. Magnetic field in the fluid frame 



(3.20) 



The magnetic field in the fiuid frame is obtained by 
substituting ([53)1 for u and ((OS)) for *F in Eq. ([XT]) . 

We get 

b = -\ \iiw + xn) - X-'e{tx,'^^,w)\ I 

(7 L L 

+ [liV - Wn) + X-'e{t X, V$, to)] X 

-(F-Tyi7)e(ex,Vf,.)}. (3.21) 

The above expression is fully general. We may 
specialize it to a perfect conductor by expressing 
e(l, X, V4>, .) via Eqs. (|XT8l) and (jXTS]) : e(^, x, .) = 
-17e(^,X, V^-,.) - iI/X)w. Using (jXTil . we get 




i(w + xn) + —df ■ V* 

An 



I v-wn- 



A2 J ■ Xn 

{V -2wn-xn^)e{tx,'^'^,-) 
I 



df ■ V* 



A 



(w + xn)w 



X 



(3.22) 



The Newtonian limit of this expression is readily obtained 
by means of Eqs. (|2.18p and (|E1[) . and after restoration 
of c^^ factors to cancel velocity terms. One gets, with 
the same notation as in Eq. (|2.46p . 



Newt.: b = 



(3.23) 

Thanks to the properties ^ • d^" = 0, x • d^* = 
[Eq. (1133)] and w ■ d'i' ^ [Eq. ([X^ ]. an immediate 
consequence of expression p.22|) is 



/:b * = 0. 



(3.24) 



Hence, like the fluid lines, the magnetic field lines are 
contained in constant 4' hypersurfaces. 
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IV. IDEAL MHD 



A. Perfect fluid model 



Prom now on, we assume that the fluid is a perfect one, 
i.e. that its energy-momentum tensor is given by 



= {e + p)u<i^u + pg, 



(4.1) 



where e is the proper energy density and p the fluid pres- 
sure. Moreover, we assume that the fluid is a simple 
fluid, that all the thermodynamic quantities depend only 
on the entropy density s and on the proper baryon num- 
ber density n. In particular, 



e — e{s, n). 



(4.2) 



The above relation is called the equation of state (EOS) 
of the fluid. The temperature T and the baryon chemical 
potential fi are then defined by 

T:=^ and ^. (4.3) 

OS on 

As a consequence of the first law of thermodynamics, p 
is a function of (s,n) entirely determined by (|4.2p : 

p = -e + Ts + fin. (4.4) 

Let us introduce the enthalpy per baryon, 



n := — fj, + 1 b, 



where S is the entropy per baryon: 

n 



(4.5) 



(4.6) 



The second equality in (|4.5p is an immediate consequence 
of dH. 



B. MHD-Euler equation 

The MHD-Euler equation stems from the conservation 
law of energy-momentum: 



V- (Tfl = 0, 



(4.7) 



where T"™ is the energy-momentum tensor of the elec- 
tromagnetic field. As is well known. 



V-T"° = F ■ j. 



On the other side, using the baryon number conservation 
V- (rtii) = 0, the term V- can be decomposed into a 
part along u, 

u-V-T^ = -nTu-dS (4.9) 
and a part orthogonal to u (see e.g. Ref. [l^ for details): 
±^,V ■ = n[u ■ d{hu) - TdS]. (4.10) 



The 2-form d{hu) is called the vorticity 2-form; 
Eq. (|4T0|) obtained first by Synge (1937) M (special 
relativity and T ~ 0), Lichnerowicz (1941) [ol (general 
relativity and T = 0) and Taub (1959) [2^ (general case) 
(see also Q) . 

From the perfect conductor relation p.3p , we have u ■ 
F ■ j = F{u,j) = 0. Hence Eq. (|4.8p has no component 
along u. Consequently, we deduce from Eqs. (j4.8|) . (|4.9p 
and (|4.10[) that the conservation law (|4.7p is equivalent 
to the system 

C^S = (4.11) 

u-d(hu)-TdS=-F-j. (4.12) 

n 

We shall call Eq. (|4TT2)) the MHD-Euler equation. 



C. Conserved quantities along the fluid lines 

For a pure rotational flow, it is a linear combination 
of the two Killing vectors [Eq. (|3.5p with w — 0] and 
every scalar field that obeys the spacetime symmetries 
is conserved along the fluid lines. This is no longer true 
for a flow with a meridional component (w 0). How- 
ever, in this case, one can derive two conservation laws 
of "Bernoulli" type, which we investigate here. 



1. Derivation 



Contracting the MHD-Euler equation (|4.12|) with the 
vector 4 leads to 



u-d{hu)-^ = F(^,j)/n, 



(4.13) 



where we have used £ • dS = since the entropy per 
baryon is supposed to respect the stationarity symmetry, 
as well as any fiuid quantity. In particular, (hu) = 
and we deduce from Cartan's identity (jB21l) that 



d{hu) ■ I = d{hu ■ ^) 



(4.14) 



Besides, from the very definition of $ [Eq. (|2.29p ]. we 
have F{^,j) = —j ■ d<I>. Using expression (|2.45a[) for 
j, along with the symmetry properties ^ • d<]? = and 
X ■ d<I> = 0, leads then to 

Fitj) = —e{t X, V/, V$). (4.15) 

fiQa 



(4.8) Thanks to Eqs. (|4l4| and (|4T5)) . Eq. (|4T3l) becomes 



Cu [hu ■ ^) 



piQan 



-e(4,x,V/,V<I'). 



(4.16) 



Since we are considering a flow which is not purely ro- 
tational, w and, from Eq. ([XT5|) . df ^ 0. Then the 
linear relation (|3.20p between d^* and df can be rewrit- 
ten as 



d^- = Cdf, 



(4.17) 
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where C is some scalar field which is necessarily a func- 
tion of /, as a consequence of the following lemma. 

Lemma 1 IJ z, p and y are three scalar fields on ^ such 
that 



dz=p(iy and dy 7^ 0, 



(4.18) 



then both z and p are functions of y, with p being the 
derivative of z with respect to y: 



and 



(4.19) 



Proof: Let us take the exterior derivative of Eq. (|4.18p 
via Eq. (jB18|) ; thanks to identities ddz = and ddy = 
[Eq. dHnJ], we get 

dpAdy = 0. 

If dp 7^ 0, this implies that the hypersurfaces of constant 
p coincide with the hypersurfaces of constant y, from 
which we deduce that p is a function of y. If dp = 0, then 
p is constant and it can still be considered as a function of 
y (constant function). Then we have dz = p(y) dy, which 
shows that z is nothing but a primitive of the function 
p{y); thus z = z{y) and p{y) = z'{y), which completes 
the proof. 

Applying Lemma [T] to Eq. (|4.17p . we obtain 



C^Cif). 



(4.20) 



Since / is preserved along the fluid lines [Eq. p.l6|) ]. we 
have of course the same property for the function C: 



£„C = 0. 

Combining Eqs. (|4.17p and p.l8|) . we get 



d$ = L> df, with D 



I 



anX 



(4.21) 



(4.22) 



From Lemma [U we have D ~ D{f) and 

C^D = 0. (4.23) 
In view of (|4.22p . Eq. (|4.16p can be rewritten as 
D 



Cu {hu ■ $,) 



-e(|,X,V/,V/). 



Now, according to Eqs. (|3.15l) and (|3.5p . 

e(^, X, V/, V/) — an w ■ dl ~ an u ■ dl — anCu I- 
Thus 

Cu {hu ■ ^) = — I = Cu 1 — ) , 



where the second equality follows from (|4.23p . We con- 
clude that the scalar quantity 



E :-- 



-hu ■ ^ 



DI 



(4.24) 



is conserved along the fluid lines 

CuE^O. (4.25) 
Thanks to Eqs. (|3.5p . we may express E as 

E^Xh{V-Wn) + —. (4.26) 
Mo 

In the limit of a vanishing electromagnetic field (/ = 
and C = 0), the conservation law ()4.25p is not hing but 
the relativistic Bernoulli theorem (see e.g. Ref. j28|). 

Repeating the same calculation, but with the Killing 
vector X instead of ^, we arrive at 



J^-a {hu ■ x) 



-e(|,X,V/,VvI') 



(4.27) 



instead of (|4.16p . Substituting Eq. (|4.17p for d* and 
making use of Eq. p.lSp . we get 



C { CI 

{hu ■ x) = — I = Cu — 
Mo V Mo 



where the second equality follows from (|4.2ip . We con- 
clude that the quantity 



L-^hu-x 

Mo 

is conserved along the fluid lines: 
CuL^O. 

Thanks to Eq. p.Sp . we may express L as 

CI 



L = \h{w + xn) 



Mo 



(4.28) 



(4.29) 



(4.30) 



The conserved quantities E and L can be considered 
as functions of /: 



E = E{f) and L 
according to the following lemma: 



L{f) 



(4.31) 



Lemma 2 // df ^ 0, any scalar field which obeys to 
the spacetime symmetries and is preserved along the fluid 
lines is a function of f . 

Proof: Let z be a scalar fleld with the above properties. 



Then C^ z = ^ ■ dz = and C^ 



% • dz = 0, which 



shows that Vz e 11-'-. Moreover, the property z = 
w ■ dz = with w (since df ^ 0) implies that 
Vz lies along the orthogonal direction to w in the plane 
n-L. The latter being generated by V/ [cf. Eq. PTT5)) ]. 
we have that dz = ad/ for some coefficient a. The 
application of Lemma [T] then completes the proof. 
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8. Newtonian limits 

To take non-relativistic limits, let us introduce the fluid 
mass density p and specific enthalpy H by 



p := nih n 



and 



H 



Jint + P 



(4.32) 



where Wb = 1.66 x 10~^^ kg is some mean baryon mass 
and fint := e — mi^n is the fluid internal energy density. 
H is related to h via Eq. (|4.5p : 



h = 7nb(l + 



(4.33) 



with _ff <C 1 at the non-relativistic limit. In view of 
(|2.18p , the expansion of Eq. (|3.6p leads to 



Newt.: A = 1 - ^grav + y, 



(4.34) 



where := w_ - w + D,^t^ sin^ 6. Substituting expressions 
((2181) . (|4:33)) and (|434l) into Eqs. (lOH)) and gJO]), we 
get the Newtonian limit of the conserved quantities E 
and L: 

E DI 

Newt.: 1 = H + $g,.av + TT + (4-35) 

mb 2 Mo™b 

r r 

Newt.: — = nr^ sin^ 61 . (4.36) 

nib p.omb 

In the absence of electromagnetic field (/ = 0), we rec- 
ognize in (|4.35p the classical Bernoulli integral. 

Besides, if we combine Eqs. (jXTT)) . ([5^ and (|iT7| . 
we recover the well known property of colinearity of the 
poloidal magnetic field and meridional velocity: 



Newt.: bp = Cn w, 



(4.37) 



L= (h+^)u-X+ — [u-{^ + ^x)]{b-x), (4.39) 
V tJ-on J Po 

where ui is defined by BO in terms of the components of 
the electromagnetic field tensor a^ 



-foi _ ^02 
^31 ^32 



(4.40) 



Using Eqs. (|E2ap - (jE2bp . we see that, within our nota- 
tions, w is the proportionality factor between the gradi- 
ents of 4> and 



(4.41) 



Combining Eqs. (|3.18p and (|4.17p . we get an expression 
of uj in terms of previously introduced quantities: 



I _ D 
CanX ^ ^C' 



(4.42) 



On Eq. (|4.4ip we see the slight shortcoming of BO expres- 
sions for E and L: if the electromagnetic field is such that 
d^* = while d<I> ^ (purely toroidal magnetic field, cf. 
Sec. I VI pp . then uj is ill defined: w — >■ oo. This corre- 
sponds to F31 = F32 = or C = [cf. Eq. (|iT7l) ]. In 
contrast, our expressions (j4.26p and (j4.30p for E and L, 
and the derivation of their constancy along the stream- 
lines, are valid even in the special case d^' = 0. Note 
however that BO formulas (|4.38p - (|4.39p give finite ex- 
pressions when w — >■ 00 ( <^=> C — > 0), since Eq. (|4.42p 
shows that 



I 



when C — >■ 0. 



INTEGRATING THE MHD-EULER 
EQUATION 



where bp is the part of b along eir) and etg\ in Eq. p.23p 



A. Explicit form of the MHD-Euler equation 



3. Comparison with BO 

The conservation laws (|4.25p and (|4.29p have been first 
established by BO 3. They have expressed E and L in 
terms of the magnetic field b in the fluid frame, but it can 
be shown that their expressions are equivalent to (j4.26p 
and (|4.30p . Actually, our derivation is slightly more gen- 
eral. Indeed, the BO expressions for E and L are0 

E^-(h+^)u-^--[u-i^ + i.x)]ib-0, (4.38) 
V Mo"- / Mo 



* Eqs. jXaSjl and l(09)l are respectively Eqs. (92) and (93) of 
Ref. @|; to show that they are equivalent to Eqs. I I4.26I I and 
14.301 1. the starting point is to use the definition 112.341 1 of / 
along with the perfect conductor expression Il3.4bl l of *F to write 



Let us first evaluate the 1-form u ■ d{hu) that appears 
in the left-hand side of the MHD-Euler equation (|4.12p . 
by means of the decomposition p. 51) of u. We first de- 
compose the 1-form u ■ d{hu) orthogonally with respect 
to the plane 11 by writing 

u-d{hu) ^ Z + aC + Px\ 

where Z is a 1-form that vanishes in 11 and the coeffi- 
cients a and /3 are determined via the properties p.37p : 
a ~ u - d{hu) ■ ^ and /3 = u ■ d{hu) ■ x- Using the Cartan 
identity (|B2ip . we get 

a = — • d{hu)] . It = — [£g {hu) —d{hu . ^)] ■ u 



= u ■ d{hu ■ ^) = w ■ d{hu ■ £). 



^ i/j is denoted —A by BO. 
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Similarly /3 = w ■ d{hu ■ x)- Hence 



Collecting ((5l| and ([5T5|) . we rewrite (jO)) as 



u ■ d{hu) = Z+[w- d{hu ■ I)] I* + • d{hu ■ x)] X* ■ 

(5.1) 

Besides, from the decomposition (|3.5p of m, we have 

u ■ d{hu) = u ■ dr + u ■ d{hw), (5-2) 

where we have introduced the 1-form 

r:=A/i(| + Ox). (5.3) 

We have, using the Cartan identity, 

M • dr = • dr + Xflx ■ dr + w ■ dr 

^X[C^r-dir-^)]+Xn[C^r-dir-x)] 



u ■ d{hu) = —A d{hu ■ ^) — Xfl d{hu ■ x) + w ■ dr 

, q 



■d/, 



(5.7) 



where we have used the property (|3.14p : e(^,x,tt,.) = 
e(l,X,i", ■) = n-^df. 

Let us employ (|5.7p to evaluate the 1-form Z acting in 
the plane 11^. Given a generic vector v e 11-'-, we have 

Z ■ V = u ■ d{hu) ■ V 

= —A V ■ d{hu ■ ^) — Xil V ■ d{hu ■ x) + dr (it), v) 

, q 



■vdf. 



(5.8) 



There remains to evaluate dr{w,v); from (|5.3|) . we have 

" " dr(w,v) ~ Xh \d^(w,v) + ridx{w,v)] , 

+w ■ dr - — 

—A d{hu ■ ^) — Xil d{hu ■ x) + w ■ dr. (5.4) where wc have used the property (dil A x)(ii', v) = 0, re- 
sulting from x-w = and x-t) — 0. By a straightforward 
calculatior0 one can show the identity 



Invoking the Cartan identity again, 

^ • d{hw) = (hw) -d{hw- 1) = 0. 
' ' 

Similarly, x ' d{hw) = 0. This shows that the 2-form 
d{h2v) acts only the 2-plane H^. By the same reasoning 
as for the 2-form G in Sec. Ill Dl we deduce that d(ft,ti;) 
must be proportional to the 2-form e(^,X: •)• 

dihw)^qe{^,X,;-)- (5.5) 
The coefficient q is determined by the Hodge duality: 



{lAxAdC){^,X,w,v) = -crdl{w,v) 
Now, the Hodge dual of relation (|2.22ap gives 

Hence 



d^{w, u) = e(|, X, v) = -^v ■ d/, 
— a an 



a \an 



where the second equality results from Eq. p.lSp . 



(5.6) 



where the second equality results from Eq. (|3.14p . Using 
the similar relation for dx(w,v), we arrive at 



dr(u;,t;) = —('^^ + n'^^)v-df. 

an 



J 



Substituting in Eq. ()5.8p . we get 

Z = -X d{hu • I) - Af7 d{hu ■ x) 

Finally, Eq. (|5.ip becomes 



Xh 



d/. 



u ■ d{hu) = [w ■ d(hu ■ |)] C + [w ■ d{hu ■ x)] X* + - 



df -Xdihu-$.)~Xnd{hu-x)- (5.9) 



Let us now evaluate the Lorentz force term on the right-hand side of the MHD-Euler equation ()4.12p . Given the 
generic form (|2.35p of i^, wc have 



F-j = ie- j) d* - (j • d$) i* + {x* ■ j) d* - (j • d^) x* + - e{t X, -J) 

a 



One may employ formula l|B5|l to express x^d^ and formula l|B3| l 
with p = 1 and g = 3 to compute ^ A (x A d^) on the quadruplet 
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Now, fromEqs. (|2.45ap and ((233)) . j-d$ = -e(|, x, V/, V'I>)/(/^ocr), and j-d* = -e(|, x, V/, V*)/(Aioo-)- Besides, 
from Eq. ^M^ . e(^, x, •, j) = -^(7^ d/. Hence 

• j = — X, V/, W) ^* + — e(C, X, W, V^/) x* + ■ j) d$ + (x* • j) - — dJ. (5.10) 



In view of Eqs. and (|5.10p . the MHD-Eulcr equation (|4.12|) becomes 



1 



-e(|,X,V/,V$) 



It) • A{hu ■ ^) — 
-X d{hu ■$,)-Xn d{hu • x) + - 



w ■ d{hu ■ x) e(^, X, V/, V^-) 

^Qan 



dJ 



df - - — - d<^> - - — - d^- - T dS* = 0. 



(5.11) 



This equation expresses the vanishing of a 1-form. The parts along ^* and x* vanish identically in the 2-plane 11-^ [cf. 
Eq. ((238)) ]. On the contrary, ah the remaining parts, being proportional to gradient of symmetric scalar fields, vanish 
identically in the 2-plane 11 = Span(^, x)- Each tangent space to ^ being the direct sum of IT and 11-'- [Eq. (|2.20p ] 
and ($*,X*) being a basis of the dual space to IT, we deduce that Eq. (|5.1ip is equivalent to the system of three 
equations: 



w ■ d{hu ■ ^) 



1 



/iocrn 
1 



-e(|,X,V/,V<i>) =0 



w ■ d{hu ■ x) e(|, X, V/, V*) = 



A dlhu ■^) + Xn d(hu • x) - - 

n 

+TdS = 0. 



£* • 7 ■ i I 

- — - d$ + - — - d^ dl 

n n jiQan 



(5.12a) 
(5.12b) 

(5.12c) 



B. Introducing the master potential 

In view of relations p.l9p and p.20p . the three linear 
forms d<I>, dvE* and df are colinear to each other. If one 
of the fields $, 4' or / is such that it gradient is non- 
vanishing, then by virtue of Lemma [T] (cf. Sec. IIV C ip . 
the two other fields can be considered as function of it. 
The standard approach in GRMHD is to privilege the 
field ^ . However this leads to degenerate equations when 
d^* = 0, which corresponds to purely toroidal magnetic 
fields or to the hydrodynamical limit (vanishing electro- 
magnetic field). The same problem occurs if one selects 
$ or / instead of ^E* (for instance selecting / leads to de- 
generate equations in the case of a pure rotational fiow). 
To be fully general, we adopt instead an approach intro- 
duced in non-relativistic MHD by Tkalich [s^ [3l| and 
Soloviev [in, namely we consider a fourth field T such 
that (i) T obeys both spacetime symmetries, (ii) dT is 
never vanishing and (iii) there exists three scalar fields 
a, /? and 7 such that 

d$ = a dT, dl* = /? dT, d/ = 7 dT. (5.13) 

The existence of T is guaranteed by the colinearity prop- 
erties (|3.19p and (|3.20p . Of course, T is far from being 
unique. The special cases mentioned above correspond to 
/3 = or 7 = 0, with dT remaining no n- vanishing. Ac- 
cording to Lemma [1] $, ^' and / are necessarily functions 



of T, with a, /3 and 7 being their derivatives: 

$ = 1>(T), vl/ = *(T), / = /(T), (5.14) 

a = $'(T), /3 = vir'(T), 7 = /'(T). (5.15) 

We call T the master potential. Using this fourth po- 
tential allows to treat all cases with finite quantities, 
whereas sticking to the three potentials and / leads 
to infinite quantities in the degenerate cases mentioned 
above. In this respect there is some analogy with the 
use of homogeneous coordinates in projective geometry: 
using only two coordinates (x, y) in the projective plane 
RP^ leads to infinite values for the "points at infinity", 
whereas adding a third coordinate, forming the so-called 
homogeneous coordinates (a;,y, z), fix this, at the price 
of some redundancy: (x, y, z) and (Az, Ay, Xz) with A 7^ 
describe the same point, as T and AT correspond to the 
same configuration. 

The master potential is conserved along any given fluid 
line. Indeed, if df ^ 0, then 7 7^ and T = u - dT = 
7-iM-d/ = by virtue of Eq. Ifd/ = 0, then m is 

a linear combination of the Killing vectors ^ and x and 
Cu T = holds according to the hypothesis (i) above. 
We conclude that in all cases 

/:„T = 0. (5.16) 

Besides, in view of (j5.13p . the perfect conductivity rela- 
tion (|3.18p is equivalent to 

a + np^^. (5.17) 
anX 
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Let us proceed by rewriting MHD-Euler system (j5.12p 
in terms of T. Thanks to Eq. p.lSp . the term 
w ■ d{hu ■ ^) in Eq. (|5.12ap can be rewritten as 

-(an)-ie(^,X,V/,V(/iM-0)- Using ([SIl, Eq. (|57[2i^ 
is thus equivalent to 



e U, X, VT, —/V{hu ■ I) + —V/ = 0. 
V Mo / 

Since 7 = 7(T) and a = Q;(T), the Leibniz rule and the 
alternate character of e allow us to write this relation as 



al 



This implies that the 1-forms dT and d{—'jhu-^+aI / hq) 
are colinear. Since dT 7^ 0, we conclude that there exists 
a scalar field a such that 



d ~^hu ■ £, 



al 
Mo 



= adT. 



Invoking again Lemma [T] we conclude that {—"fhu ■ ^ + 
a//Mo) must be a function of T, I](T) say. Expressing 
M • i via Eqs. ([331) and (|2T3l) - (f2J6l) . we get 

0.1 dl 
S(T) = -7/iu, • I + — = -i\h{V - Wn) + — . (5.18) 
Mo Mo 
I 



Applying a similar argument to the second equation of 
the MHD-Euler system ()5.12|) leads to the existence of a 
function A(T) such that 



A(T) = -/hu-x - 



Mo 



-fXh{W + XQ) - 



Mo ' 



(5.19) 



As for any function of T, the quantities S and A are 
conserved along any given fluid line, in consequence of 

(EUl). 

Note that if d/ ^ 0, one may choose T — f, leading 
to the following values [cf. Eqs. (|4T7l) . (g^l), (14261) and 



T = f 




E 

A: 



E 
L. 



(5.20) 



Hence, for this choice of T, E and A are simply the 
Bernoulli-like quantities E and L introduced by BO Q 
and discussed in See. IIV CI 

If d^ 7^ 0, the choice T = is allowed, leading to 




E = E/C 
A = L/C, 



(5.21) 



where uj = -D/C [ef. Eq. (gUl)] 



C. The master transfield equation 



Having shown that the first two equations of the MHD-Euler system (|5.12p leads to the conserved quantities E and 
A, let us focus on the third equation, namely Eq. (|5.12e|) . Taking account of ()5.13|) . we can rewrite it as 



\d{hu-Cl + >^^A{hu-x) + 



dT 



Mocn 



■dH-TdS' = 0. (5.22) 



Differentiating expressions (jS.lSp and (|5.19l) yields (the prime stands for the first derivative of a function of T) 

/ a'I\ a 

dE = E' dT = -7d(/iM ■$,)-{ -y'hu ■ $, dT + — d/ 

V MO / MO 

dA = A' dT = jd{hu ■ x) + ( I'hu ■ x - — I dT + ^ d/, 

V Mo / Mo 



from which we get 
\d[hu ■i) + Xn d{hu ■x) = - 



nA' - E' + — (a' + Qp') + -/'xh{v - 2wn - xn^) 

Mo 



dT 



a + ni3 
Mo 



dl 



To treat the entropy term T dS in Eq. (|5.22|) . let us assume that S* is a function of T: 

S = S{T). 



(5.23) 



(5.24) 



Actually (|5.24p is mandatory if there is a non-vanishing meridional flow. Indeed in this ease df ^ and since S is 
conserved along the fluid lines [property (|4.1ip ]. Lemma [2] of See. IIVC H is applicable and gives S = S{f), i.e., via 
(|5.14p . 5* = S'(T). If df ~ (pure rotational motion), then we may consider that (|5.24p is a supplementary hypothesis 
in our framework, set to integrate the MHD-Euler equation. Note that a homentropic fluid (5 = const, throughout 
the fluid) satisfles 
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Substituting Eqs. (|5.23p into Eq. (|5.22[) . we notice that terms in dl cancel each other thanks to the relation (|5.17l) . 
Using (|5.24p to set dS = S' dT, we are then left with 



nA' - s' + —{a + np') + i\h{v - 2wn - xn'^) 

Mo 



TS' MT = 0. 



(5.25) 



Since by hypothesis dT ^ 0, this equation is equivalent to the vanishing of the term in braces. Let us express all the 
pieces in term of T. Writing d<I> = adT and d* = /3dT [Eq. ([STT^ ] in Eqs. (j2.45b|) and (|2.45cp . we get 

aC ■ 3 + PX* ■ 3 = ^—\{Vli^ + 2WaP-Xa^)A*T + {P^dV + 2al3dW~a^dX) ■ VT 

/iQCr L 

+ [V/Sfi' + W{a'/3 + a/3') - Xaa'] dT • VT - - [{Wfi - Xa)'^. + {Wa + VfSy^^] |,(5.26) 

a J 

where A* is the operator that generalizes p.47p to the relativistic case: 



A*T aV„ ( ^V^T 



Besides, setting d/ — ^dT [Eq. (|5.13p ] in Eq. (|5.6p . we have 

h'~f 



-A*T + 7d 



VT + -7'dT- VT 



(5.27) 



(5.28) 



Let us substitute Eqs. (|5.26p and (|5.28p into the term in braces in Eq. (|5.25p and express its vanishing. We get, after 
multiplication by an? /h, 



AA*T + - 
h 



77 



72d ( - ) - — {p^dV + 2a(3dW - a^dX) 



VT 



— [VP13' + W{a'p + aP') - Xaa'] \dT-VT 



an 



QA' - s' + —{a + nfi') + j'xh{v - 2wn - xn^ 

Mo 
In 

-7A71 {'tf^ + fl'f^) + [{WI3 - XaYf^ + {Wa + V^Y^^] = 



TS' 



where 



A 



{V/3^ + 2WaP - Xa^). 



(5.29) 



(5.30) 



We shall call Eq. (|5.29p the master transfield equation. The qualifier transfield stems from the fact that it corresponds 
to the component of the MHD-Euler equation along dT, which is transverse to the magnetic field: b ■ dT = 0, as it is 
easily verified on the expression p.22p of b, taking into account (j5.14p and (j5.16p . Note that the term in A/7 in the 
third line of Eq. (|5.29p is regular, even when 7 = 0, as we shall show in Sec. IVI C II Besides, note that for circular 
spacetimes (e.g. Kerr spacetime), the last line of Eq. (|5.29p vanishes identically [cf. Eq. (|2.2ip ]. The master transfield 
equation was first written in the Newtonian case by Soloviev (1967) [U, as we shall discuss in Sec. IVI Al 

Given the metric [hence the covariant derivative operator V, the coefficients V, W, V and a, and the twist scalars 
and and the six functions a{T), /3(T), 7(T), S(T), A(T), and S{T), the master transfield equation (|5.29p 
constitutes a (non-linear) second-order PDE for T. Indeed, all the remaining quantities (/, A, il, n, h, T) that appear 
in Eq. (|5.29p . although not functions of T, can be computed once T is known, as we will show. 

First of all, by combining Eqs. (|5.18p and (|5.19p . we get 



VA~WT. = jaXhn {VI3 + Wa) 

XJ: + WA^ -faXh + —{Xa - W^/3). 



(5.31) 
(5.32) 
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Combining these two equations and using (j5.17p to express a + in terms of /, we get 



Tl 

I = —[{Xa-Wl3)J: + {Wa + V(5)A]. (5.33) 



Extracting \h from Eq. (|5.32p and substituting Eq. (|5.33p for / leads to 

1 



Xh 



(TjA 



j^X^ + WA) - ^(aA + /3S) 
Mo" 



(5.34) 



Extracting from Eq. (|5.3ip and substituting the above values of / and Xh, we obtain the expression of the fluid 
angular velocities in terms of conserved quantities and h/n: 



Ha^-f^iVA - WJ:) + aaiPT. + aA) 
/io^7'(^S + WA) - a/3(/?E + aA) ' 



Besides, from the relation (|3.15p . we have 



1 ^ 72 ^ 

w w^ — ^d/- V/= ^dT- VT, (5.36) 



so that the 4-velocity normalization relation (j3.6p can be written as 



^2 

1 + -^rdT -VT = X^(V -2Wn- Xn"^). (5.37) 



Substituting Eq. (|5.35p for ft and writing A = {Xh)/h with Xh given by (|5.34p . we get, after some rearrangements, 



(^^ + ;b dT . VTj - ^(XE^ + 2W^A VA') + — ^^-^(/^E + aA)' = 0. (5.38) 

By means of the identity 

cr(/3I] + aA)2 = [{Xa - WP)T. + (VP + Wa)Af + [VP^ + 2Wal3 ~ Xa^){XY.'^ + 2WY.A - VA^), 
which follows solely from a = XV + , Eq. (j5.38p can be recast in the alternative forrrQ 

(^<T + ^dT- Vrj - ^ (XE^ + 2M^EA - VA^) + ^t^^^ [(^" ' ^^)^ + (^^ + Wa)Af = 0. (5.39) 

This equation is called the poloidal wind equation. Given the metric factors V, W , X and cr, the functions a(T), /3(T), 
7(T), E(T), A(T) and ^(T), expression dOU)) for A, as well as the EOS h = h{s,n) with s = S{T)n [cf. Eq. (ITO ]. 
the poloidal wind equation can be solved to compute n once T is known. Then, from n we get h via the EOS and A 
via Eq. (|5.30p . Once n, h and A are known, we can compute / via Eq. (|5.33p and D, via Eq. (j5.35p . The meridional 
velocity w is obtained via Eq. p.lSp with df = 7 dT and the velocity coefficient A via Eq. p.6p . 

An equivalent point of view is to consider that the fundamental equations to be solved are Eqs. (|5.29p and (|5.39p 
which constitute a coupled PDE system for the two unknowns (T,ri). Indeed, given the metric, the EOS and the 
six functions a(T), /3(T), 7(T), E(T), A(T) and S{T), solving this system provides a solution of the MHD-Euler 
equation and Maxwell equations, the electromagnetic field tensor F and electric 4-current j being deduced from T 
via Eqs. ([OS)) . ^A5^, ([Fl3)) and ((03| . 

VI. SUBCASES OF THE MASTER TRANSFIELD EQUATION 

The master transfield equation (|5.29p . coupled with the poloidal wind equation (|5.39p . describes the most general 
MHD equilibria in generic (noncircular) stationary and axisymmetric spacetimes. We shall now specialize it to various 
cases and make the link with results obtained previously in the literature. 



7 A combination -S + HA of Eqs. JESJ and JSHJ may be used ^^^-^^ g^l and JSJOj from u u = -l. 
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Newtonian limit 



At the Newtonian limit, as given by Eqs. ((2?T8| . ([23T|) . and (|434)) . the expressions (fSTTS)) and ((5?T9| of the 

streamline-conserved quantities I](T) and A(T) reduce to 



E = 7?«b ( 1 + i? + $grav 



V 

~2 



A ~ jnihr sin" 0fi — 



Mo 



whereas the master transfield equation (|5.29p reduces to 

„2 



AA*T -—dn 

n 



2 • 2 

-r sm ( 



VT 

2 ri 



77 



^/3' dT • VT 



71 

mb 

with the expression (|2.47|) for the operator A* an41 



ilA' - S' + — (a + nf3') + 7'mb 
Mo 



A 



n/S'^ 
Mo^b' 



r5' ^ = 0, 



(6.1) 
(6.2) 



(6.3) 



(6.4) 



Performing the proper changes of notatioro, one can check that Eq. (|6.3p coincides with the first equation in the 
system (II) of Soloviev [s^ . After Soloviev work, the Newtonian transfield equation has been re-obtained by many 
authors for the special case in which T = ^I* [cf. (|5.21|) ] (e.g. [l^, [ill, [13, [lOl, [U; cf. [s^ for an extended 
discussion and [53[ for a recent study). The equation is then known as the generalized Grad-Shafranov equation (see 
Sec. EE] below). 

The Newtonian limits of expressions (j5.33p and (|5.35p for / and il are 



T ^ I 2-2 

7 = — [ ajr sm ' 



1 



7A 



sin^ I 



PA 
mb 



— ap 
Mo 



(6.5) 
(6.6) 



To get the latter expression, we have used the Newtonian limit /3S] -I- aA = mb/37, resulting from (|6.ip - (|6.2p . Equa- 
tions (|6.5p and (|6.6p coincides with respectively the second and first equations in Eq. (1.18) of Soloviev article [32| . 

Finally to get the Newtonian limit of the poloidal wind equation (|5.39p . we rewrite the term h^a — ArS^/7^ by 
means of the Newtonian expressions (|2.18p of a and X, along with ([4.33^ : 



7^ 



2 • 2 
r sm I 



m^(l + 2i7)- (1 



J]2 

2*I*grav) 2 



2m1r^ sin^ 9 



mb7 



where the last equality results from 1 — (I]/mb7)^ = (1 + I]/mb7)(l — S/™b7) — 2(1 — I]/mb7)- Accordingly, the 
Newtonian limit of the poloidal wind equation (|5.39p is 



^dT • VT + 2r2 sin^ 6{h^ $„.av + 1 



— ) 



A 

mb7 



A + 72 / /3A 



^on^b A^ V™b7 



ar'^ sir? ( 



0. (6.7) 



This equation is not exhibited in Soloviev work [33 |. It can however be recovered by combining Soloviev Eqs. (1.5) 
and (1.25) and expressing as ^I'^r^ sir? 9 + w_- w with substituted by (|6.6p and w-why (|5.36p . In the special 
case where T = ^E", one can check that Eq. (|6.7p coincides with Eq. (14) of Heyvaerts and Norman [5l[ (called the 
Bernoulli equation by these authors). 

I 



* Note that in taking the Newtonian Hmit of 1)5. 30p . the term 
Xa^ = sin^ 9a^ is relativistic and therefore must be disre- 



garded. 
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B. Pure rotational flow 

The case of a pure rotational flow corresponds to 

w = ^ d/ = ^ 7 = 0. (6.8) 
Then Eq. ([5T7)) yields 

(6.9) 

whereas Eqs. (|5.18p and (|5.19p reduce to 
al . . /?/ 



and A 



Mo 



(6.10) 



If a 7^ (i.e. d$ 7^ 0) or /? ^ (i.e. d^- ^ 0), Eqs. (lOl) 
and (|6.10|) imply that ft and / are functions of T (for a, 
/3, S and A are all functions of T): 



n = and / = /(t). 



(6.11) 



Taking into account Eq. (|6.9p and 7 = 0, the expression 
(|5.3Gp for A becomes 



3.12) 



To express the master transfield equation (|5.29|) in the 
case 7 = 0, we shall first evaluate the term which is 
divided by 7 in Eq. (|5.29p . namely 

A := nA'-j:' + —{a' + np') + j'xh{v-2wn~xn^). 

Mo 

(6.13) 

To this aim, we shall first suppose 7 7^ and, in a second 
stage, take the limit 7 — > 0. Since / = ^(T) when 7 = 6 
[Eq. (|6.1ip ]. we may write 



/ = /o(T)+7a, 



(6.14) 



where /o(T) is a function of T and a describes the behav- 
ior of / as 7 — )■ 0. For instance, if /3 7^ (i.e. d^* 7^ 0), 
explicit values of Iq and a are deduced from Eq. (|5.19p : 

/o(T) = -Mo||^ and a ^ ^io^{W + XQ). 



Substituting expression (I6.14p for / into Eqs. (|5.18p and 
([ETO]) . we get 

S = —lo + iG and A = -— /q + 7H, (6.15) 
Mo Mo 



^ The link between Soloviev notations |32|| and ours is r r sin 0, 
s mbA/n, mb7, "I'' <-> f)/^/Jm, lo Hr^ sin^ 

/ o I/^/Jms, S •<-> S/m^, A A, _B a/^/JTo, and C7 •(-!> 
{■mb7)-i(S + QA//3)-l. 

with 

G := A/i(V^-P^l]) + — and iJ := A/i(M^ + Xm - ^. 

Mo Mo 

(6.16) 

From Eq. (|6.15p and the fact that S, A, a, /3, 7 and /q 
are all functions of T, it is clear that G = G(T) and 
H — H{T). Then, using successively Eqs. (|6.15p . (|6.16p 
and (|5.17p . we may write expression (|6.13p as 



A = -i{nH' -G' + — 
Mo 



crnX 



w+«y) 



(6.17) 

Taking the limit 7 0, we have D, = n{T) [Eq. ((6TTt ] 
and a = -rj/3 [Eq. (US])], so that a' + flf3' = and 
ri H'-G' = -K'-n'H with K{T) := G~nH. According 
to (jOBI and dnH), we have K = Xh{V - 2Wn - ATJl^). 
Now for 7 = 0, expression (|3.6p for A reduces to 



A= (F-2W^17-Xl]2)-i/2^ 



so that 



A'(T) = - = h^/v -2wn- xn^. 

A 



(6.18) 



(6.19) 



Since 



ilH' -G' + —{a' + nf3') = -K' - n'H - — f7'/3 
Mo Mo 



-K' - xh{w + xn)n' , 



the limit 7 ^ of Eq. ((6l7| is 
A 



lim 

7— >o 7 



-K' - Xh{W + XVl)Vl' 



II' 
MocnA 



(6.20) 



Thanks to the relations (|6T2l) . (|6T8l) . (|O0|) and (|4?5|) . the master transfield equation (|5?29l) becomes 



/3 



+M0O- Ue +p) 



^ A*T + ^ d ( ) • VT + /?2(VK + ATf])!]' dT • VT + /<{/'- ^ [(M^ + Xny<f^ + {V - Wil)"^^ 



/3 



x^{w + xn)n' 



nTS' 



0. 



3.21) 
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Let us now take the pure rotational limit of the wind equation, in the form (|5.38p . From Eq. (|6.15p . /3S + aA = 
7(/3G' + aH), so that, in view of (|6.9|) . 

lim -{/3T. + aA) = /3{G - QH) = ^K. 
7— >0 'y 

This property, along with (|6.12[) . implies that for 7 = the wind equation (|5.38|) reduces to 



V - 2wn - xvfi 

which is nothing but the square of Eq. (|6.19|) . Consequently, the wind equation is trivially satisfied in this case. 

In conclusion, for a pure rotational fiow, one should prescribe five functions of the master potential: /?(T), il(T), 
/(T), K{T) and S'(T) and solve for the transfield equation (|6.2ip for T. In that equation, the matter quantities e, 
p, n and T are given by the EOS from the knowledge of S and /i, the latter being deduced from r2(T) and i^(T) via 
Eq. (|6.19p . We shall discuss further the pure rotational fiow below (Sec. IVIC H and IVID2)) . 



C. Expression in terms of ^P: generalized Grad-Shafranov equation 



Let us assume that 7^ 0. We may then choose T = '5 as the primary variable [cf. (|5.2ip ]. This is actually the 
choice performed by most (all ?) of previous relativistic studies, disregarding the case = (toroidal magnetic 
field or hydrodynamical limit, to be discussed in Sec. I VI Dl and IvTeI . Let us first consider the pure rotational fiow, 
in order to make the link with the original Grad-Shafranov equation. 



Pure rotational flow 



For a pure rotational flow with T = 4*, Eq. (|6.1ip become 



n = ni'ii) and / = /(*) 



3.22) 



In the Newtonian regime, the property 57 = 51(^') is known as Ferraro's law of isorotation [33, while the result 
/ — /(5') has been obtained by Chandrasekhar (1956) [13] • The transfleld equation ()6.2ip becomes [cf. ()5.2ip and 

dHHl)] 



{V - 2wn - xn^)A*^ + d{v - 2wn - xn^) • + (vk + xn)n' d* • v* + / 



I' - 



w + xn 



+^j.a(7 \{e+p) 



K' {xn + w)n' 



V - 2wn - xn^ 



nTS' 



0. 



(6.23) 



This PDE has to be solved for ^P, once the four functions 5^(5'), ^(^'), A'('I') and S{'^) are prescribed. The enthalpy 
field h is given by Eq. (|6.19p which remains unchanged. 

Equation (|6.23p is a relativistic generalization of the so-called Grad-Shafranov equation [35l - l38| (sec also Chap. 16 
of the recent textbook [l^). At the Newtonian limit [cf. expressions (|2.18p and (|4.32p ] and in coordinates (t, r, 0, </?) 
of spherical type, Eq. (|6.23p reduces to 



Newt.: 



A** + //' + Mo?-^ sin^ e [p {K'/rrih + fin'r^ sin^ 6) - nTS'] = 0, 



whereas Eq. (pl^ reduces to [cf. and 



Newt.: H + $grav ^^r"^ sin^ 



rrih 



1. 



3.24) 



3.25) 



Let us recall that the Newtonian expression of A* is given by Eq. (|2.47p . Modulo the change from spherical to 
cylindrical coordinates, Eqs. (|6.24p and (|6.25p coincide with respectively Eqs. (3.3) and (3.2) of Maschke and Perrin 
[40|. The limiting case / = [pu re p oloidal magnetic field, cf. Eq. (|3.23p ]. p = const, = const and S = const has 
heeii treated by Ferraro in 1954 [4l[. It has been extended to / ^ and 51 7^ const, still maintaining p = const, by 
Chandrasekhar in 1956 [s^] (using the function P := /{r'^ sin^ 9) instead of ^). Plasma physicists Grad and Rubin 
[35} and Shafranov [s^] have considered in 1958 the non-rotating limit (fi = 0) of Eq. (|6.24p (see Chap. 16 of (39|). 
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Coming back to the relativistic case, the special case 7 = 0, = const and S = const has been discussed by 



Bonazzola et ah [2J| and Bocquet et al. [25[. Note that contrary to what is claimed in Ref. [2^, has not to be 
a constant: it can be any function of 4' [Eq. ()6.22|) ]. Most relativistic studies have focused on the case it) 7^ (flow 
with a meridional component) and barely discussed the limit w = presented above. In particular, it is claimed in 
Ref. ^3 that if k; — ^ 0, the magnetic field b cannot have a toroidal component (i.e. / = 0). We see no support of this 
since any choice seems to be allowed for the function 1(5') in the equations presented above. 



2. Generic flow 



For a generic flow (i.e. with some meridional component), we have, according to (|5.2ip . a 
7 = C^^, with C 7^ since 7^ 0. The expression ()5.30|) for A becomes then 



A = 



1 

CP 



where M is the poloidal Alfven Mach-number 



1 - 



V - 2Wlu - Xuj' 



M2 



-w, /3 = 1 and 



(6.26) 



This name is justified by the Newtonian limit [cf. (|4.33p ] 



Newt.: 



VA,] 



VA.: 



(6.27) 



(6.28) 



^A,p is the poloidal Alfven velocity, bp being the poloidal magnetic field [cf. Eq. (|4.37p ]. The expression of M as the 
ratio of the norm of w to VA.p justifies the name poloidal Alfven Mach-number given to M . 

Setting Y. = E/C and h = L/C [cf. ([OT]) ]. the transfield equation ([Og]) speciahzed to T = * is 



V - 2Wlo - Xu? 



M2 



M2 



{W + Xuj) 



c 



1 

/iQCrn 



{dV - 2ujdW - w^dX) 



V* 



A 



VlL' -E' + — {C'{n -co)- CJ) 



TS' 



I 



XnC + ^ + {{W + XuY^^ + (V^ - Wijj)'€^ = 0. 



3.29) 



In this equation, all the primes denote derivatives with respect to Equation (|6.29p is called the generalized 
Grad-Shafranov equation, since it can be considered as an extension of the Grad-Sliafranov equation (|6.23p to the 
case of a non-vanishing meridional fiow. The generalized Grad-Shafranov equation has been derived for Minkowski 
spacetimes by Camenzind (1987) [53 and Heyvaerts & Norman (2003) [55|. It has been extended to weak gravitational 
fields by Lovelace et al. (1986) j50| and to the Schwarzschild spacetime by Mobarry and Lovelace (1986) [l^. The 
case of the Kerr spacetime has been first considered by Nitta et al. (1991) [ll| for pressurclcss matter and Bcskin 
and Pariev (1993) [l^l for non-vanishing pressure (see also [l|). Finally the general case of noncircular stationary 
axisymmetric spacetimes has been treated by loka and Sasaki (2003) [ll,[l^. Note however that they have not written 
the generalized Grad-Shafranov equation explicitly as Eq. (|6.29p . but have kept the ^* ■ j, X* ' 3 and q terms, as in 
Eq. (|5.25p . They have replaced these terms, leading to A*^, only when taking the Newtonian limit. In particular, 
it is not apparent in their work that the Grad-Shafranov equation is singular at the Alfven surface, where the term 
in factor of A*\E' vanishes (see below). Besides, as stated in the Introduction, their approach appeals to a (2-|-l)+l 
foliation of spacetime, whereas ours does not require any extra structure. 
Regarding the poloidal wind equation (j5.39p . it reads for = T, 



1 



d-ii ■Vi' + ah^ - XE^ - 2WEL + VL^ 



2AP - 1 [{V- Woj)L -{W + Xuj)E\ 



(Af2 - 1)2 



V - 2Wuj - Xw2 



0, 



wher( 



m 



Af2 := 



M2 



V-2Wu:-X(^'^ C^{V -2Wuj - Xui^) n' 



(6.30) 



(6.31) 



Note that M is not necessarily positive, contrary to M 
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Equations (j6.29p and (|6.30p form a system of two equations for the two unknowns n), once the five functions Ld{^), 
C(5'), E{'i>), L(^) and are prescribed, as weh as the EOS. In these equations, / is expressed via Eq. (|5.33p : 



and n via Eq. (|5.35p recast as 



/xo (V - Wuj)L -{W + Xlo)E 
'C M^-V^ + 2Wuj + Xw2 



M^{VL - WE) - auj{E - ljL) 
M^{XE + WL) - a{E - wL) ' 



(6.32) 



(6.33) 



As a check, one may verify that Eq. (16.321) coincides with Eq. (24) obtained by Beskin & Pariev [12| for the Kerr 
spacetimj"! and that Eq. (|6.33|) coincidcjj with Eq. (110) obtained by Camenzind [i^ for the Minkowski spacetime. 
It can also be recovered for general circular spacetimes by combining Eqs. (38), (39) and (41) of Ref. [Zsj . 

The generalized Grad-Shafranov equation (|6.29p is singular for = V— 2Wuj — Xlo"^, or equivalently, for A'P = 1. 
This condition defines the so-called Alfven surface (see e.g. Refs. [ll. I43l - l46j for an extended discussion). The term 
— 1 also appears at the denominator in the the poloidal wind equation (|6.30p or in expression (|6.32p for /, 
but this does not make these equations singular at the Alfven surface, thanks to the simultaneous vanishing of the 
corresponding numerator [l| . 



D. Toroidal magnetic field (d* = 0) 

The complementary case of that treated in the previous subsection is 

d^- = <=^ /3 = 0, 



(6.34) 



where the equivalence follows from the very definition of /3 given in Eq. (|5.13p . Then, the expression p.22p for the 
magnetic field in the fluid frame reduces to 



a 



w • w \ 1 

iw + xn)$+ {V -wn-=^j x- jiw + xn)w 



(6.35) 



Strictly speaking, this field is not purely toroidal, except at the Newtonian limit or when w = 0. By a slight abuse of 
language, we shall however refer to the case d\l/ = as the toroidal magnetic field case. 



1. Generic case 



With /3 — 0, the perfect conductivity relation (|5.17p reduces to 



Consequently, the expression (|5.30p of A becomes 
The master transfield equation (|5.29p reduces to 



XP 



AA*T + Y 



P 



-dX 



fioaXn ) ri \l 



liQcrXh 



(-A<^4 + W'S'^) = 0, 



(6.36) 



(6.37) 



Mo 



(6.38) 



11 The link between notations of Ref. [11 and ours is as follows: 12 ^j^^ ^^j^ between notations of Ref. M and ours is a o a; and 
O a/X, uj o -W/X, « aX^jX, ro^ o X, / -7/2, ij2 o <t 

* «• Ato>I'/2, r] «• « £ «• EjC and L « L/C. 
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whereas the poloidal wind equation (|5.39p becomes 

+ ^ dT • VT^ - \ (XS^ + 2H^SA - VA^) + —(2+ , ) = 0. (6.39) 



2. Pure rotational flow 

In the particular case of a pure rotational flow (7 = 0), 
the transfield equation (j6.38p reduces to [cf. Eq. (j6.2ip 
with /3 = 0] 



//'+/iocr \{e + p) 



K 



-nTS' 



0, 



(6.40) 

with A'(T) obeying to Eq. (|6.19p . In the present case, 
a /? = 7 = 0, i.e. d$ = 0, d* = and d/ = 0. If the 
electromagnetic field is not vanishing, a natural choice 
for T is 



T = /. 



Then /' = 1 and Eq. ([QUI reduces to 



(6.41) 



1{ 



x^{w + xn)n' 



nTS' 



3.42) 

Given the functions ri(/), K{I) and S'(/), this equation 
has to be solved in /. Note that this is not a PDE in / 
and that the matter quantities n, e, p and T are to be 
computed via the EOS from S and h, the former being 
deduced from K{I) and via Eq. (pl^ : 



/v(/) 



^y- 2iyji(/) - xi}{iy 



3.43) 



In the special case ^{I) = const and S{I) = const, we 
recover equations obtained by Kiuchi and Yoshida (2008) 
pT^ El- At the Newtonian limit, the special case — 
const, S{I) = const and K'(I)/K{I) = const has been 
contemplated by Miketinac (1973) |42| . 

It is worth underlining that in the case considered here, 
i.e. a pure rotational flow and a pure toroidal mag- 
netic field, (i) the spacetime has to be circular (provided 
that the fluid and the electromagnetic field are the only 
sources in the Einstein equation) [2^ [23| and (ii) the 
twist functions and "^-^j whose vanishing is equivalent 
to circularity, do not appear in Eqs. (|6.42p - (|6.43p . 



The link between notations of Ref. [27 
92 CT, \/ 92/91^12 -n- I, u -n- anh - 
-U.oK'{I)/K{I). 



and ours is g\ — g/cr, 
: o-(e +p), 4:1tK{u)/u ^ 



E. Hydrodynamical limit 

1. Generic case 

The hydrodynamical limit (no electromagnetic field) is 
easily taken by setting a = (i.e. d$ = 0), /3 = (i.e. 
d\I/ = 0) and / = in the equations obtained so far. 
In particular, the streamline-conserved quantities (jS.lSp 
and (|5.19p reduce to 

E = 7 £; with E = \h{V - Wn), (6.44) 
K^-iL with L = \h{W + Xn), (6.45) 

where the second equalities in each line follow from 
Eqs. (|4:26| and (jiTSOl with / = Bes i des, E q. (lOOl re- 
duces to A = 7^ and thanks to (i04|) - (|05l) . Eq. (|03l) 
reduces to ^ = ^{VLL' — E'). Accordingly the master 
transfield equation (|5.29p becomes 



72A*T + 72-d 
h 



VT + 77'dT ■ VT 
A(fiL' - E') + TS'] - jXn ("^^ + il'tfy^) 



0. 

(6.46) 

On its side, the poloidal wind equation ()5.39p reduces to 
J_ dT • VT + ah^ - XE^ - 2WEL + VL^ = 0. (6.47) 

S. Flow with meridional component 

If the meridional fluid velocity is not vanishing, d/ 7^ 
and a natural choice for the master potential is T = /. 
Then 7=1 and Eqs. (|6.46p and (|6.47p become 



A * ,. n , f h 
A*/ + -rd - 
h V 71 



[\{^IL' - E')+TS'] 



-An ("4 ^U^) = 



(6.48) 



■ d/ • V/ + ah^ - XE^ - 2WEL + VL^ = 0.(6.49) 



These equations are to be supplemented by (i) Eq. 
expressing A in terms of / and fl [via Eq. (|5.36p ] and (ii) 
the hydrodynamical limit of Eq. (|5.35|) . which reads 



VL-WE 
XE + WL' 



(6.50) 



It is then clear that, given the metric, the three func- 
tions £:(/), L{f) and S{f) and the EOS h = h{S,n), 
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T = T{S,n), Eqs. (|08)) - (|09)) form a system of cou- 
pled PDE for {f,n). Solving this system provides a 
solution of the Euler equation for a rotating flow with 
meridional component. In the case of circular spacetimes 
='^x = 0)iEq. ((08l) has been written first by An- 
derson (1989) [Hil and Bcskin & Pariev (1993) [13, with 
n extracted from Eq. (|6.49p and substituted in (|6.48l) , so 
that the system reduces to a single equation for /. An 
equivalent formulation has been developed recently by 
Birkl et al. [53| for the barotropic case {S = 0), using 
the function tp := — J E{f) df instead of /. 

In the Newtonian limit [cf. Eqs. (|6.3p and (|6.7p ]. the 
system becomes 

1 



A*/--dn- V/ 



n 



T T ' 

— -r^ sin2 e {E' - TS') 

?71b 



(6.51) 



4j d/ • V/ + 2r2 sin2 e{H + $g,av + 1 - — 
L 

mb 



3.52) 



Note that we have substituted Vt by the Newtonian limit 
of n = L/{mhr^ sin^ 6). Equation (|S3T|) is called 

Stokes equation. For an incompressible fluid, it coincides 
with Eq. (7.5.11) of Batchelor treatise [H^. For a com- 
pressible fluid, we recover Eq. (1.107) of Beskin textbook 
[l| or Eq. (15) of Eriguchi et al. [53 • Consequently, we 
shall call Eq. (j6.48p the relativistic Stokes equation. 



3. Pure rotational flow 



For a pure rotational flow, 7 = and Eqs. (|6.46p and 
([07)) reduce to 



flL' - + ^5' = 
A 

ah^ = XE'^ + 2WEL - VL^. 
Let us assume that 

n = n{T). 



(6.53) 
(6.54) 

(6.55) 



In Sec. IVIBI we have seen that this condition is manda- 
tory if a ^ or /3 ^ 0. In the absence of electromagnetic 
field, ()6.55p is fulfilled for rigid rotation (f2 = const) and 
for non-rigid one, it may be imposed by a proper choice 
of T, for instance T = fl. It is then natural to intro- 
ducers in Sec. lyiB] K = K {T) = E - flL = h/X [cf. 
Eq. (|6.19p ]. so that Eq. (|6.53p can be written as 



K' + m' - ^S" = 0. 

A 



(6.56) 



A wide class of equilibrium configurations is obtained by 
assuming 



T 
A 



T(T), 



(6.57) 



with T is a function of T such that 

f = -f-Vt'. 
K 



(6.58) 



Thanks to Eqs. (|6?57|) - (|638l) . Eq. (|636l) can be written 
as 



(In^)' + —n' = 0. 
A 



where 



^ = /2(T) -.^K-TS 



h — TS jjL 
A " A' 



(6.59) 



(6.60) 



jjL being the baryon chemical potential introduced in 
Eq. (|4.3p . the last equality resulting from Eq. (|4.5p . 

If f2 is constant (rigid rotation), Eq. (|6.59l) yields the 
first integral 



/i = const. 



(6.61) 



Note that in this case, Eq. (|6.58p implies the so-called 
relativistic isothermal condition : 



T = const. 



3.62) 



If n' ^ 0, Lemma n of Sec. |IVC1J implies that L/K 
must be a function of i7, J^(ri), say. Since L is expressible 
as (|05)) and K ^ h/X with the value (pl^ for A, we 
have 



w + xn 



V - 2w^i - xn^ 



(6.63) 



Given the function J^{^), the above equation has to be 
solved in r2. At the Newtonian limit, it leads to a solution 
of the form = ri(r sinO), i.e. satisfying to the Poincare- 
Wavre property With L/K = J^{n), Eq. ((^3^)) is 

integrated to 



In fi+ / J^(f^) dfl = const. 



and Eq. to 



const, 



(6.64) 



(6.65) 



generalizing the isothermal condition (|6.62p to the case 
of differential rotation. In the case T = (or S = const), 
we recognize in Eqs. (|6.6ip and (|6.64p the standard first 
integrals governing rotating relativistic stars (see e.g. [6l| 
or J8j and references therein). For the finite temperature 
case, we recover results of Ref. 1631. 



VII. SUMMARY AND CONCLUSION 

We have formulated GRMHD for stationary and ax- 
isymmetric spacetimes in the most general case, i.e. non 
assuming circularity (as in Kerr spacetime). Moreover, 
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wc have based our approach on geometric quantities de- 
fined solely in terms of the spacetimes symmetries (rep- 
resented by the two Killing vectors ^ and x), without 
relying on any coordinate system or any extra structure 
(such as a (2-|-l)-|-l foliation). This provides some insight 
on previously introduced quantities and leads to the for- 
mulation of very general laws, recovering previous ones as 
subcases and obtaining new ones in some specific limits. 
To our knowledge, the new results obtained here are: 

• the expression (|2.35p of the electromagnetic field 
tensor F entirely in terms of the two Killing vector 
fields and three scalar fields, independently of any 
coordinate system; 

• the derivation of the conservation laws for the 
Bernoulli-type quantities E and i in a covariant 
manner and in the most general case, including that 
of a purely toroidal magnetic field disregarded in 
the original BO derivation Q; 

• the fully covariant master transfield equation (j5.29|) 
governing the most general MHD equilibria in 
generic (i.e. noncircular) spacetimes, generalizing 
Soloviev non-relativistic equation (3^ : 

• the explicit form (|6.29|) of the covariant Grad- 
Shafranov equation for noncircular spacetimes; 

• the equation (|6.38|) governing MHD equilibria with 
purely toroidal magnetic field in stationary and ax- 
isymmctric spacetimes; 

• the relativistic Stokes equation (|6.48|) governing hy- 
drodynamical equilibria of flows with meridional 
components in stationary and axisymmetric space- 
times. 

The relativistic master transfield equation (|5.29p is prob- 
ably the most important outcome of the present study. 
Beyond the aesthetic feature of having a single equation 
governing all MHD equilibria, reducing to the relativistic 
Grad-Shafranov and Stokes equations in certain limits, 
the value of this equation resides in its potentiality to 
lead to solutions that cannot be obtained by merely set- 
ting T = or T = /, as already shown in the Newtonian 
regime [63j . 

In this article, we have focused on the derivation of 
the equations governing MHD equilibria and of conser- 
vation laws. In order to solve the obtained equations, 
there remains to choose the streamline-conserved func- 
tions a(T), /3(T), 7(T), I](T), A(T) and 5(T) and to 
specify some boundary conditions on them. An example 
of numerical resolution of the Grad-Shafranov equation 
(case T = is provided by Ref. [l^. Finding station- 
ary and axisymmetric GRMHD solutions provides initial 
data for dynamical stability studies of magnetized neu- 
tron stars (see e.g. [63 - l67j ). 

As a final remark, let us point out that we have hardly 
used the axisymmetric character of the Killing vector x 



(i.e. the fact that it is a generator of a S0(2) group ac- 
tion), so that most results presented here would remain 
valid for any other type of spatial symmetry, like for in- 
stance translational symmetry. 
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Appendix A: Lie derivative 

The Lie derivative is the natural operator to express 
symmetries under a 1-parameter group action. It mea- 
sures the change of a tensor field along the orbits of the 
group action. More precisely, any regular vector field v 
can be regarded as the generator of a 1-parameter group 
action on Then, given a local coordinate system 

(x") adapted to v, i.e. such that v" — (1,0,0,0), the 
Lie derivative along v oi a tensor field T of type (fc, i) 
is the tensor field of the same type, whose components 
are the derivatives of T's components with respect to the 
parameter associated with v (i.e. the coordinate x'^): 

(^«T^r-"V...ft = ^or"^-%,...^,. (Ai) 

In an arbitrary coordinate system, this formula becomes 

i=l 
I 

+ E^"-"ft..-.../J.^ft«^ (A2) 

i=l i 

In particular, for a scalar field /, 

C.f^ v^d^f, (A3) 

for a vector field w, 

(£„ wr = v^'d^.w'^ - w^d^v-^, (A4) 
for a 1-form a;, 

(£„ U})a = vf^dfj^UJa + UJf^daV'^, (A5) 
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and for a bilinear form T (such as the metric tensor g or 
the electromagnetic field -F), 

TUi = v^d^^Ta^n + T^pdc^v^ + Ta^^dpv''. (A6) 

From formula (|A4[) . note that the Lie derivative of w 
along V is nothing but the commutator of the vector fields 
V and w : 

CvW = [v, w]. (A7) 

Note also that in formulas (|A2[) - (jA6p . one may replace 
the partial derivative operator d by the covariant deriva- 
tive operator V associated with the metric g. This stems 
from the symmetry property of the Christoffel symbols. 
In particular, Eq. (jA6p can be written 

r)„^ = t.^V^T„^ + T^pV^v^' + T^^V^w^. (A8) 



Appendix B: Differential forms and exterior calculus 

Given a integer p satisfying < P < 4, a p-form is 
a tensor of type (0,p) which is fully antisymmetric. By 
convention, a 0-form is a scalar and a 1-form is a linear 
form. A differential form of rank p is a field of p-forms 
over Differential forms play a special role in the 

theory of integration on a manifold. Indeed the primary 
definition of an integral over a manifold of dimension p 
is the integral of a p-form. At a given point x € 
the set of all 1-forms is the dual vector space 



to the vector space tangent to ^ at x. More 

generally the set A^{^) of all p- forms at x is a vector 
space of dimension (^). In particular, the dimension of 
the space of 4-forms is 1: all 4- forms are proportional to 
each other. 

We assume that the manifold ^ is orientable, i.e. 
that there exists a continuous, nowhere vanishing, 4-form 
field. We may then introduce the Levi-Civita alternating 
tensor e (also called metric volume element) as the dif- 
ferential form of rank 4 such that for any vector basis 
(cq,) that is orthonormal with respect to g, 

6(60,61,62,63) = ±1. (Bl) 

If ^ is orientable, there are actually two such 4-form 
fields, opposite to each other: picking one of them is 
making a choice of orientation on Having a -I- sign 
(resp. — sign) in Eq. (jBl[) defines then a right-handed 
basis (resp. a left-handed basis). The components of e in 
a given right-handed basis (not necessarily orthonormal) 
are 

eap-yS = V^[a, P, 1, S\ , (B2) 

where g is the determinant of the components (gap) of 
the metric tensor in the considered basis and [a, /3, 7, 5] 
stands for 1 (resp. —1) if (a,/?, 7, (5) is an even (resp. 
odd) permutation of (0, 1, 2, 3), and otherwise. 

Two algebraic operations are defined on differential 
forms: the exterior product and Hodge star. The exte- 
rior product associates to any p-form A and any g-form 
B the {p -I- q)-form A A B defined by 



AAB{vi,...,Vp+q) := — (-l)'=('"U(t;^(i),...,v<,(p)) X B(v<,(p+i),...,v<,(p+g)), (B3) 

where Vi,...,Vp+q are generic p -\- q vectors, &p+q is the group of permutations of p -\- q elements, (~1)'=('^) is the 
signature of permutation a and and x denotes the multiplication in R. In particular, if A and B are 1-forms, the 
exterior product is expressible in terms of tensor products as 

AAB^A(E)B-B(E)A (1-forms). (B4) 

If A is a 1-form and B is a 2-form, then 

AAB{vi,V2,V3) = {A ■ vi) B{v2,V3) + {A ■ V2) B{v3,vi) + {A ■ V3) B{vi,V2). (B5) 



The Hodge star operator relies on the Levi-Civita ten- called the Hodge dual of u), and defined by 
sor e: it associates to every p-form w, a (4 — p)-form ★a;, 

0- form : {*u)af3yS =■ ^ f-afi-iS (B6) 

1- form : {i^u>) afi-, = e^„^^ (B7) 
1 

2' 

1 

—L 

6 

1 

24" 



2- form : {-ku})ajj = -uj^^ e'"'^^ (B8) 

3- form : (*u;)„ = Itu^^p (B9) 

4- form : ^u; = T^w^^p^ e'^'P" . (BIO) 
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Notice that, for any p-form, 

★ ★w=(-l)f+iw (Bll) 
and that, for any couple {a, h) of 1-forms. 

* (a A b) = e(a, b, ., .) and ★ [e(a, b, ., .)] = a A b, 

(B12) 

where a (resp. b) is the vector associated to a (resp. b) 
by the metric [cf. Eq. (|^ ]. 

Being tensor fields, the differential forms arc subject to 
the covariant derivative V and to the Lie derivative £„ 
discussed above. But, in addition, they arc subject to a 
third type of derivation, called exterior derivation. The 
exterior derivative of a p-form field a; is a (p + l)-form 
field denoted do;. In terms of components with respect 
to a given coordinate system (a:"), Auj is defined by 

0- form : (du;)^ = d^uj (B13) 

1- form : [du})^^ — daOJp — dpuja (B14) 

2- form : (da;)^^^ = daUip^ + dpuj^a + d^^ap (B15) 

3- form : {du))apjs = daUJpjs - dpui^sa 

+djU!SaP - dsUlaP-y- (B16) 

It can be easily checked that these formulas, although ex- 
pressed in terms of partial derivatives of components in a 
coordinate system, do define tensor fields. Notice that for 
a scalar field (0-form), the exterior derivative is nothing 
but the gradient 1-form. Notice also that the definition 
of the exterior derivative appeals only to the manifold 
structure. It does not depend upon the metric tensor g, 
nor upon any other extra structure on Besides, as 
for the Lie derivative expressions (jA2[) - (|A6[) . all partial 
derivatives in formulas (|B13p - (|B16p can be replaced by 
covariant derivatives V thanks to the symmetry of the 
Christoffel symbols. 

A fundamental property of the exterior derivation is to 
be nilpotent: 

ddw = 0. (B17) 

A p-form u} is said to be closed iff da; = 0, and exact iff 
there exists a (p — l)-form cr such that a; = dcr. From 
property (jB17p . any exact p-form is closed. The Poincare 
lemma states that the converse is true, at least locally. 

With respect to the exterior product, the exterior 
derivation obeys to a modified Leibniz rule: if a is a 
p-form and b a (/-form, 

d(aAb) ==daAb+(-l)PaAdb. (B18) 

If p is even, we recover the standard Leibniz rule. 

The exterior derivative enters in the well known Stokes ' 
theorem: if 2? is a submanifold of of dimension d e 
{1,2,3,4} and has a boundary (denoted dT)), then for 
any {d — l)-form u;, 

* u}^ I du3. (B19) 

JdV JV 



Note that dT> is a manifold of dimension d — 1 and da; 
is a d-form, so that each side of (|B19p is a well defined 
quantity, as the integral of a p-form over a p-dimensional 
manifold. 

A standard identity relates the divergence of a vector 
field V to the exterior derivative of the 3-form v ■ e (see 
e.g. Appendix B of Ref. [2^): 

d(t;-e) (V-^;)e. (B20) 

Another very useful formula where the exterior deriva- 
tive enters is Cartan identity, which states that the Lie 
derivative of a p-form a; (p > 1) along a vector field v is 
expressible as 

C^u) = v-du: + d{v-u)). (B21) 

Notice that for a 1-form, Eq. (jB2ip is readily obtained 
by combining Eqs. (|A5p and (|B14p . 

Appendix C: Relation between the scalar fields <1> 
and 'i/ and the electromagnetic 4-potential 

The electromagnetic 4-potential A is not an observable 
and may not necessarily obey the symmetries (|2.28p of 
the electromagnetic field F = dA. However, for each 
Killing vector, one can find a gauge transformation such 
that the 4-potential obeys the corresponding symmetry. 
We demonstrate this for stationarity and axisymmctry. 

Stationarity implies that C^F = dA = dC^ A = 
0. If ^ is simply connected, then the Poincare lemma 
implies that there exists a single- valued scalar fj. such that 

A = dn, but this quantity will be nonzero in general. 
However, there exists a class of gauge transformations 
A' = A + diy such that £^ A' = C^{A + diy) d{fi + 
C^v) ~ provided that the scalar v satisfies ^ + v = 
const. This differential equation may be integrated along 
the integral curves of the timelike Killing vector ^. This 
procedure eliminates the time-dependent part of A. 

Axisymmctry implies the relation C-^^ F ~ C-^ dA' = 
dC^ A' ~ which, by virtue of the Poincare lemma, im- 
plies the existence of a single-valued scalar fi' such that 
C-^A' — d/i', but this quantity will again be nonzero 
in general. However, there exists another class of time- 
independent gauge transformations A" = A' + dv' such 
that A" = {A' + diy') = d(//' + v') = with 
a time-independent scalar v' (obeying C^^v' = Q by as- 
sumption) that satisfies the equation /j,' -I- v' ~ const. 
This differential equation can be integrated along the in- 
tegral curves of the axial Killing vector The resulting 
gauge transformation eliminates the non-axisymmetric 
part of A' while maintaining its stationarity. 

This permits one to work in a gauge class within which 
the 4-potcntial A is stationary and axisymmetric. From 
the Cartan identity, ^ • dA -I- d(^ - A) = £| A = 0, one 
then has ^ ■ F ^ — d^t and similarly x ' F ~ —dA^p. 
Comparing these two equations to (|2.29p and (j2.30p al- 
lows one to identify At with $ and A^p with ^P, up to some 
additive constant, thereby demonstrating Eq. (|2.3ip . 
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Appendix D: Kerr-Newman electromagnetic field Xil = x* • M. Using p.39p and (|2.13|) - (l2.16p . we get 



The Kerr-Newman solution describes a charged rotat- 
ing black hole. In Boyer-Lindquist coordinates {t, r, 6, ip), 
its electromagnetic field is |68| 



4:n{r'^ + cos 



2a\2 



{P hdt + Rhde) , (Dl) 



where Q is the total electric charge, a J/M the 
reduced angular momentum of the black hole, P := 
(r-^ — cos^ 6) dr ~ a'^r sin 26 d9 and R := a(a^ cos^ 9 — 
r^)sin^6'dr -|- ar{r^ + a^) sin26'd0. Since the Kerr- 
Newman spacetime is circular, ^* = dt and x* = 
[cf. Eq. <^Z3l] . The comparison with ([OS)) leads to 



$ = - 



/ = 0. 



^J■oQ r 
Att 4- cos^ 9 ' 



Mo I 



ar am 



Att + a? cos^ 9 ' 
(D2) 



At the non-rotating limit (a = 0), this reduces to 
Reissner-Nordstrom solution: 



47r r ' 



* = 0, 



/ = 0. 



(D3) 



Appendix E: Component expressions with respect to 
a coordinate system 

In coordinates (x") = {t,x^ ,ip) adapted to sta- 
tionarity and axisymmetry (cf. Sec. Ill Bp . the compo- 
nents of the Killing vectors are ^" = (1,0,0,0) and 
= (0,0,0,1) [cf. Eq. (pm)) ]. From Eq. the 
2-form e(4, •) is expressible as 



1. General spacetimes 



(El) 



The components Fap of the electromagnetic field are 
given by F Fapdx"" ® dx^ = F\af}\dx°' A dx^, 
where \ol(3\ means that the summation is limited to 
a < p. Accordingly, substituting Eqs. p.40p and (|ET|) 
into Eq. (|2.35p leads to 



Foa 
Fl2 





1 r 



a e {1,2} 
a e {1,2} 



9i$(-^6 + 1^X2) - d2<i>{-X^i 
di^{W^2 + Vx2) - d2-^{Wii 



(E2a) 
(E2b) 
(E2c) 

Wxi) 
-Vxi) 
(E2d) 



Regarding the fluid 4- velocity m, the coefficients A and 
f2 in the decomposition p.Sp are given by A = ^* • it and 



a 



(E3) 
(E4) 



From and the fact that = and = (a e 

{1,2}), we get 

w"" = u"". (E5) 

From (|3.5p and (|E3P - (|E4p . the remaining components of 
w are 



1 



= -{Xia~WXaK 

a 



(E6) 
(E7) 



From the property w E 11-'-, we get wq = and 

""^3 = w^x^ = 0> hence the covariant components of w: 



Wa = (0, Ua- Xi^a+^Xa), 0) 



(E8) 



The covariant components of u arc given by uq — UfiS,'^ 
and U3 = u^x'^j so that we may write 

i-X{V -Wn), Ua, XiW + Xil)) . (E9) 



2. Circular spacetimes 

In the circular case, we may choose coordinates 
(t, x^,x'^, (p) so that the surfaces orthogonal to 11 arc the 
surfaces {t = const, ip — const}. Then (|2.23p holds and 
we have = .ga^C^ = .9a0 = and Xa = gaf^x'' = .9a3 = 
0: 



L=Xa^ 0. 



(ElO) 



Accordingly, the components (jE2p of the electromagnetic 
field simplify to 



/ 
































(Ell) 



and relations (IE3I)-(IE4I) reduce to 



^ X and u' 
whereas (|E6P - (|E7p become 

w° = w^ = 0. 



xn. 



(E12) 



(E13) 
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